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ABSTRACT 

The  Ocean  Thermal  Energy  Conversion  (OTEC)  uses  the  low 
thermal  energy  potential  available  from  ocean  temperature 
gradients.   A  method  is  presented  to  analyze  such  systems  and, 
for  this  purpose,  a  comprehensive  simulation  is  developed. 
The  simulation  includes  parasitic  power  requirements,  losses 
due  to  interconnecting  lines,  and  heat  exchanger  pressure 
drops.   Cost  functions  are  included  and  numerical  optimization 
is  employed  to  obtain  optimal  designs  based  upon  minimum  cost. 
The  analysis  is  converted  to  a  computer  code  and  coupled  to 
the  COPES/CONMIN  optimization  code  to  facilitate  a  fully- 
automated  design  where  the  computer  makes  the  design  decisions 
and  performance  trade-off  studies.   The  final  product  is  an 
optimum  power  system  module  design  for  the  designated  net 
electrical  output  required  and  the  specified  system  and  design 
constraints . 

Preliminary  results  are  presented  for  a  range  of  system 
power  levels.   Optimum  designs  are  obtained  and  compared  for 
systems  in  which  either  titanium  or  aluminum  tubes  are  used 
in  the  heat  exchangers. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Ocean  Thermal  Energy  Conversion  (OTEC)  is  a  concept 
using  the  low  thermal  energy  potential  available  from  the 
ocean  temperature  gradient  that  exists  between  warm  surface 
Ocean  water  and  cold  water  in  deep  ocean  regions. 

The  idea  of  converting  the  stored  ocean  energy  to  useful 
power  originated  with  French  physicist  Jacques  d'Arsonval  in 
1881  [Ref .  1] .   It  was  nearly  a  half-century  later  that  the 
technical  feasibility  of  ocean  thermal  energy  conversion 
could  be  demonstrated.   In  1926,  George  Claude  used  an  open 
cycle  power  system  to  extract  heat  from  surface  water  for 
indirect  conversion  of  the  thermal  energy  of  a  working  fluid. 
Operating  at  a  low  pressure  the  working  fluid  was  used  to 
drive  a  turbine  providing  electrical  power  generation. 

Though  Claude's  limited  power  system  produced  only 
22  kilowatts  of  electricity  while  requiring  approximately 
80  kilowatts  of  power  to  drive  its  equipment,  it  stirred 
the  scientific  and  research  community  to  consider  the 
attractiveness  of  ocean  thermal  energy  conversion  [Ref.  2]. 

Claude  called  for  immediate  action  on  his  ocean  thermal 
power  system,  because  of  the  Federal  Oil  Conservation  Board's 
dire  predictions  that  the  United  States  had  only  six  years 
of  oil  production  remaining.   Obviously  the  dire  predictions 
ascribed  to  by  the  Federal  Oil  Conservation  Board  did  not 
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come  true,  but  the  oil  crisis  of  that  period  heightened 
scientific  interest  in  extracting  energy  from  the  ocean. 

Now,  55  years  later,  the  United  States  is  faced  with  an 
energy  crisis  because  of  increasing  industrial  and  social 
dependence  on  foreign  petroleum.   Dwindling  supplies  and 
erratic  price  hikes  have  rekindled  interest  in  ocean  thermal 
energy  conversion,  since  it  utilizes  an  inexhaustible  supply 
of  fuel. 

Currently,  the  United  States  Department  of  Energy  is 
attempting  to  develop  the  necessary  technology  and  demon- 
strate the  feasibility  of  large-scale  OTEC  power  systems. 
However,  there  are  major  engineering  development  problems 
which  must  be  solved  before  OTEC  can  be  standardized  and 
become  a  viable  source  of  electrical  power  generation. 

The  single  controlling  factor  which  creates  troublesome 
technical  encounters  is  low  thermal  power  system  efficiency 
(one  to  four  percent  depending  upon  parasitic  power  require- 
ments) .   Because  the  heat  energy  used  by  OTEC  must  be  extracted 
from  a  small  ocean  temperature  difference,  extremely  large 
volumes  of  surface  water  must  pass  through  a  proportionately 
sized  evaporator  to  provide  sufficient  indirect  heat  energy 
to  convert  the  working  fluid  into  vapor  to  drive  a  turbine - 
generator  for  electrical  power  generation.   Concurrently,  to 
convert  the  turbine  exhaust  to  a  saturated  liquid,  completing 
the  closed  cycle,  a  condenser  having  compatible  heat  absorp- 
tion capacity  must  be  employed. 
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Economic  handling  of  the  volume  of  fluids  required  for  the 
heat  absorption,  expansion,  and  heat  rejection  phases  of  the 
cycle  requires  close  scrutiny  of  evaporator,  turbine,  condenser, 
and  pump  design  to  minimize  the  parasitic  losses  with  respect 
to  the  generated  electrical  output.   Because  of  the  low 
thermal  efficiency,  relative  to  nuclear  or  fossil  fuel-fired 
power  plants,  the  margins  for  design  and  operating  error  in 
OTEC  plants  will  be  narrow. 

With  the  advent  of  high-speed  computers,  numerical 
methods  for  solving  these  complex  engineering  problems  with 
multiple  design  variables  and  constraints  are  now  possible. 
The  case  for  utilizing  an  optimizing  scheme  for  not  just 
one  system  component,  but  rather  the  complete  power  generation 
cycle,  can  easily  be  made.   In  effect,  it  would  serve  as  a 
systems  analysis  tool,  to  optimize  component  design  and  cost, 
relative  to  a  specific  electrical  output  or  to  enable  compari- 
son and  evaluation  of  competing  OTEC  designs. 

B.   OBJECTIVES 

The  objectives  of  this  work  are  to  develop  a  computer 
code  for  the  Ocean  Thermal  Energy  Conversion  (OTEC)  power 
system  and  to  couple  the  analysis  to  a  numerical  optimization 
code  to  provide  an  optimum  system  design  capability,  considering 
both  performance  and  economics. 

This  would  create  an  optimum  modular  design  relative  to 
a  specified  objective  function  for  a  desired  net  electrical 
output,  such  as  a  25  MW  (net)  power  system.   Such  a  design 
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would  permit  construction  of  higher  capacity  power  systems 
using  the  optimized  modules  as  substations  of  the  total  power 
plant.   Cost  savings,  improved  plant  performance,  redundancy, 
and  reliability  could  be  the  immediate  beneficiaries  of  such 
a  venture . 

C.   OVERVIEW  OF  THE  OTEC  POWER  SYSTEM  ANALYSIS 

To  analyze  the  closed-cycle  OTEC  power  system,  the 

fundamental  relationships  of  heat  transfer,  fluid  mechanics 

and  thermodynamics  are  used  to  simulate  a  variety  of  system 

component  designs,  which  form  the  basis  of  the  power  system 

algorithm.   The  scope  of  this  analysis  will  be  limited  to 

the  OTEC  power  system  and  sea  water  systems  only.   Mooring 

systems,  power  delivery,  hull,  and  cold  pipe  design  will  not 

be  addressed. 

The  performance  analysis  will  be  divided  into  four 

sequential  sections  as  shown  in  Figure  1,  and  discussed 

in  detail  in  subsequent  chapters  of  this  thesis. 

Input  parameters  (design  constants)  for  the  power  cycle 

analysis  will  include: 

.  Required  net  electrical  output. 

.  Salt  water  inlet  temperature  to  the  evaporator 
and  condenser. 

.  Length  of  hot  and  cold  salt  water  pipes. 

.  Heat  exchanger  tubing  material  (aluminum  or  titanium) 

.  Heat  exchanger  tube  orientation  and  profile. 

.  Pump  mechanical  and  motor  efficiencies. 
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Turbine  mechanical  efficiency. 

Generator  mechanical  and  electrical  efficiency. 

Biofouling  control  factor. 

Piping  absolute  roughness. 

Projected  annual  inflation  rate  for  aluminum  heat 
exchanger  retubing. 
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Figure  1.   Power  System  Sequential  Analysis 
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II.   POWER  CYCLE  DESCRIPTIONS 

A.  INTRODUCTION 

This  chapter  will  provide  a  brief  description  of  the  OTEC 
power  system.   First,  looking  at  the  ideal  Rankine  cycle, 
the  fundamental  thermodynamic  concepts  will  be  enumerated. 
Then  the  deviations  from  the  ideal  cycle  will  be  presented, 
creating  the  configuration  assumed  for  the  present  cycle 
analysis  which  will  be  amplified  in  detail  by  follow-on 
chapters . 

B.  IDEAL  OTEC  RANKINE  CYCLE 

The  closed-cycle  OTEC  concept  is  based  upon  a  Rankine 
power  cycle  that  is  driven  by  the  low  thermal  energy  potential 
available  from  the  ocean  temperature  gradient  that  exists 
between  warm  surface  water  and  cold  deep  water  in  ocean 
regions.   The  power  cycle  consists  of  a  working  fluid  circu- 
lation pump,  evaporator  (heat  absorption),  turbine  (expansion), 
and  condenser  (heat  rejection),  as  shown  in  Figure  2.   The 
majority  of  current  OTEC  designs  are  based  upon  ammonia  as 
the  working  fluid  --  a  design  decision  that  is  adopted  for 
this  analysis. 

Figure  2  also  illustrates  an  ideal  OTEC  Rankine  cycle, 
plotted  on  temperature -entropy  coordinates.   In  the  ideal 
cycle,  the  low  pressure  working  fluid  (state  point  1)  is 
isentropically  pumped  to  the  evaporator  operating  pressure 
(state  point  2) .   The  working  fluid  (ammonia)  is  then 
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converted  to  a  saturated  vapor  in  the  evaporator  by  indirect 
heat  energy  exchange  from  warm  surface  ocean  water  (state 
point  3) .   Mechanical  power  is  generated  by  isentropic  expan- 
sion of  the  saturated  ammonia  vapor  through  the  turbine 
Cstate  point  4) . 

After  exiting  the  turbine,  the  wet,  low-pressure  vapor  is 
converted  to  a  saturated  liquid  in  the  condenser  by  indirect 
heat  absorption  from  cold  ocean  water  (state  point  1) , 
returning  the  cycle  back  to  the  working  fluid  circulation 
pump . 

C.   ACTUAL  OTEC  RANKINE  CYCLE 

In  actuality  there  are  numerous  deviations  from  the  ideal 
cycle  which  must  be  considered  in  this  analysis.   These  are: 

(1)  Turbine,  generator  and  pump  efficiencies. 

(2)  Pressure  drops  in  evaporator  and  condenser  (tube- 
side  and  shellside) . 

(3)  Pressure  drop  across  moisture  separator. 

(4)  Elevation  change  and  frictional  losses  in 
piping:   (a)  re-flux  pump  piping,  (b)  piping 
from  circulation  pump  to  evaporator. 

(5)  Evaporator  outlet  quality  (85  to  95%) . 

(6)  Moisture  separator  outlet  quality  (99  to  99.51). 
The  deviations  from  the  ideal  Rankine  cycle  described 

above  are  depicted  in  the  flow  diagram  and  temperature-entropy 
plot  of  Figure  3.   In  the  actual  OTEC  Rankine  cycle,  the  low 
pressure  working  fluid  (state  point  1)  is  pumped  up  to  the 
evaporator  operating  pressure  by  the  ammonia  circulation  pump 
with  an  adiabatic  efficiency  (state  point  2) .   The  working 
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fluid  (ammonia)  is  then  converted  to  a  wet  vapor  with  an 
evaporator  outlet  quality  (85-95%)  acting  under  a  shellside 
pressure  drop  (state  point  3) .   Evaporator  outlet  vapor  then 
passes  through  a  moisture  separator  to  improve  vapor  quality 
(99-99.51)  creating  a  pressure  drop  (state  point  4). 
Mechanical  power  is  generated  by  the  expansion  of  the  moisture 
separator  outlet  vapor  through  the  turbine  with  an  adiabatic 
efficiency  (state  point  5).   After  exiting  the  turbine,  the 
wet,  low  pressure  vapor  is  converted  to  a  saturated  liquid 
in  the  condenser  acting  under  a  shellside  pressure  drop 
(state  point  1) ,  returning  the  cycle  to  the  working  fluid 
circulation  pump. 

This  figure  forms  the  thermodynamic  basis  for  the  OTEC 
power  system  analysis  which  follows. 
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Ill .   EVAPORATOR  AND  MOISTURE  SEPARATOR 

A.   INTRODUCTION 

Several  heat  exchanger  concepts  have  been  proposed  for 

closed-cycle  OTEC  systems.   Among  these  designs  are: 

.  Conventional  shell  and  tube  heat  exchanger. 

.  Plate  type  heat  exchanger. 

Within  these  basic  concepts,  variations  in  design  have 

been  proposed,  including: 

.  Orientation  of  tubes  (horizontal  or  vertical). 

.  Heat  exchanger  tube  material  (i.e.,  titanium, 
aluminum)  . 

.  Method  of  tube  enhancement  (i.e.,  fluted,  porous 
coatings) . 

.  Location  of  tube  enhancement  (i.e.,  internal  and/ 
or  external) . 

.  Location  of  the  vapor  separator  (i.e.,  internal 
or  external) . 

.  Location  of  the  heat  exchangers  relative  to  the 
sea  surface. 

.  Method  of  biofouling  control. 

The  analysis  to  be  presented  for  the  evaporative  heat 

exchanger  will  be  based  on  the  following  design  characteristics 

.  Single-pass  shell  and  tube  heat  exchanger. 

.  Internal  vapor  separator  with  a  gravity  drain  to 
evaporator  inlet. 

.  Horizontal  orientation  of  tubes  with  an  equi- 
lateral triangle  or  square  tube  profile. 

.  Smooth  plain-tube  configuration  (no  enhancements). 
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.  Tube  material  (titanium  or  aluminum  based  on  a 
3Q-year  life-cycle  criterion) . 

.  Biofouling  control  based  upon  an  achievable 
fouling  factor. 

.  Heat  exchanger  centerline  located  on  sea  surface. 

As  an  overview  of  the  evaporator-moisture  separator 

analysis,  the  following  major  steps  in  the  algorithm  are 

listed  in  order  of  their  execution  (numbers  in  parentheses 

refer  to  equations  developed  in  the  subsequent  analysis)  : 

.  Specification  of  system  constants  (see  I.C.). 

.  Initialization  of  design  variables  (D.V.). 

.  Tube  length. 

.  SW  velocity  through  hot  pipe. 

.  Inner  diameter  of  hot  pipe. 

.  Tube  outer  diameter. 

.  SW  velocity  through  evaporator  tubes. 

.  Inner  diameter  of  NH3  piping. 

.  Inner  diameter  of  NH3  re-flux  piping. 

.  Tube  profile  pitch  ratio. 

.  Salt  water  mass  flow  rate  (1) . 

.  Total  number  of  tubes  (2) . 

.  Total  heat  transfer  surface  area  (3) . 

.  Assume  an  initial  salt  water  bulk  temperature  (6) ,  and 
ammonia  heat  transfer  coefficient  (9) . 

.  Overall  heat  transfer  coefficient  (4) . 

.  Number  of  transfer  units  (11). 

.  Heat  exchanger  effectiveness  (13) . 

.  Salt  water  outlet  temperature  (15) . 
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.  Revised  bulk  temperature  (16)  ;  iterate  with  (6) . 

.  Amount  of  heat  absorption  (17) . 

.  Log  mean  temperature  difference  (18)  . 

.  Film  temperature  (19) . 

.  Initial  ammonia  mass  flow  rate  (21)  without  the 
effects  of  moisture  separator. 

.  Initially  assume  state  point  1  thermodynamic 
properties  are  ideal  (21)  . 

.  Thermodynamic  pump  work  (23) . 

.  Tube  profile,  flow  parameters  across  the  tube 
bank  (24,  etc .) . 

.  Tube  sheet  diameter  (30)  . 

.  Evaporator  shellside  pressure  drop  for  two  phase 
flow  (33). 

.  Moisture  separator  pressure  drop  (38) . 

.  Properties  at  state  points  3  and  4  (39-41). 

.  Revised  ammonia  mass  flow  rate  and  velocity  (50) 
includes  the  effects  of  the  moisture  separator; 
iterate  with  (31)  . 

.  Revised  ammonia  heat  transfer  coefficient  (51,  etc.); 
iterate  with  (9) . 

.  Heat  exchanger  cost  analysis. 

In  the  following  section,  the  basic  steps  summarized 

above  will  be  described  in  detail. 

B.   ANALYSIS  OF  THE  EVAPORATOR  AND  MOISTURE  SEPARATOR 
1 .   Salt  Water  Mass  Flow  rate,  KQsw 

The  salt  water  mass  flow  rate  through  the  hot  pipe 
must  be  equivalent  to  the  flow  rate  through  the  evaporator 
(assuming  no  leakage) 
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r 


'yy\    I  Hot  ?}?£.)-   YY\  z  (E  M  Po  £/\  ro  *) 


and 


rn  =/i,<  a  v 


(i) 


where 


A   =  cross -sectional  area  of  the  hot  pipe. 

V  =  salt  water  velocity  through  hot  pipe. 

,/^vv  =  density  of  salt  water  evaluated  for  an 

average  hot  pipe  salt  water  temperature. 

As  previously  stated,  the  diameter  of  the  hot 
pipe  and  salt  water  velocity  are  among  the  initializing 
conditions  of  the  optimization  and  will  be  treated  as  design 
variables . 

2.   Total  Number  of  Evaporator  Tubes,  fsjt 


Using  equation  (1) ,  it  follows  that 


m1  «/iw  ffdi   \A  Nv 


(2) 


-4 


where  /&sn  ■  salt  water  density  evaluated  at  the  average 

bulk  temperature  initially  assumed  as  the  hot 
pipe  salt  water  temperature. 
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dl       =  tube  inner  diameter. 
N-t   =  the  number  of  tubes  required  to  maintain 

the  mass  flow  rate  for  an  average  salt  water 
velocity  per  tube. 

The  total  number  of  tubes  can  be  determined  by  solving 
Eq.  (2)  for  Mt- 

The  diameter  of  the  tube  and  average  salt  water 
velocity  per  tube  are  initialized  for  the  analysis  and  will 
be  treated  as  design  variables  by  the  otpimization  code. 

3.   Total  Evaporator  Heat  Transfer  Surface  Area  (Outer),  A< 

Having  determined  the  number  of  evaporator  tubes,  the 
total  heat  transfer  surface  area  can  be  determined  using 
initializing  values  of  outer  tube  diameter  and  tube  length. 

For  tubes  without  extended  surfaces 


A±  =  l^dc  Lt  Nt  (3) 


As  previously,  the  outer  tube  diameter  and  tube  length 
are  initializing  conditions  and  will  be  treated  as  design 
variables . 

4.   Overall  Heat  Transfer  Coefficient,  Li. 


The  quantity  "U"  provides  a  measure  of  the  total 
thermal  resistance  in  the  flow  path,  based  on  either  inside 
or  outside  surface  area. 

This  analysis  will  be  based  on  the  value  of  U  for  the 
outside  surface  area  derived  from  Eq.  (3) . 
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Using  a  resistance  analysis,  assuming  one  dimensional 
(radial)  heat  flow, 


■sw 


Tsw   A  A   Tfsw   a  a    Twl 


w2 


vVv-^V\A^-^WV-^x^V-i^yVV 


Tnh- 


the  overall  heat  transfer  coefficient  may  be  expressed  as 


u.- 


where 


( 


Ac       Acn      do  In  do  J  di 


R* 


i 


(4) 


tiH3 


hsw  =  tubeside  heat  transfer  coefficient. 
Rfsw  =  salt  water  fouling  heat  transfer  resistance. 

K   =  thermal  conductivity  of  the  tube  material. 
d0)di    =  outer  and  inner  tube  diameter. 
RfWrt  =  ammonia  fouling  heat  transfer  resistance 

(assumed  to  be  negligible)  . 
77*,^L'=  ou'ter  and  inner  total  fin  efficiency  (for 
plain  tube  analysis,  total  fin  efficiency 
equals  1)  . 
/\Q     =  total  outer  surface  area  (including  fin 
and  bare  tube) . 
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Ai      =  total  inner  surface  area  (including  fin 
and  hare  tube) . l 


where     Af?i;     =  total  inner  fin  surface  area. 

Ai       ■  total  inner  surface  area  (including  fin 
and  bare  tube) . 
Afri^     =  total  outer  fin  surface  area. 
/^   =  total  outer  surface  area  (including  fin 

and  bare  tube) . 
}[,■      =    fin  efficiency  of  single  interanl  fin. 
7V   =  fin  efficiency  of  single  external  fin. 


a.   Tubeside  Reynolds  Number,  R<ej 

Since  the  heat  transfer  coefficient  correlations 
for  the  evaporator  and  condenser  are  dependent  on  tubeside 
flow,  Reynolds  number  must  be  calculated. 

The  tube  Reynolds  number  is  defined  as 


p     ,/2s»   VW  dl  (5) 


*Note  that  this  analysis  will  hereafter  consider  smooth 
plain  tube  configurations  only. 
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where  /ttsv*   -   dynamic  viscosity  of  salt  water. 
fsv*    =  density  of  salt  water. 

Initially,  properties  are  evaluated  for 

Tbulk   =  7s*j  (MlcTJ  (6) 

Reynolds  numbers  greater  than  2300  will  be  indica- 
tive of  turbulent  flow  [Ref.  3].   Transition  flow  was  considered 
laminar  for  numerical  evaluation. 

b.   Salt  Water  Heat  Transfer  Coefficient,  n^^ 

The  simple  empirical  relation  proposed  by  Sieder 
and  Tate  [Ref.  3],  expressed  as 

was  used  for  laminar  heat  transfer  in  tubes  as  defined  by 
Eq.  (5). 

Nusselt  and  Prandtl  numbers,  htUji   and  \r  ,  are 


defined  as 


K 


S(Ai 


where  zf^,  £p       and  Ks.u  (dynamic  viscosity,  specific  heat, 
and  thermal  conductivity)  of  salt  water  are  evaluated  at  salt 
water  bulk  temperature. 
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The  effect  of  the  viscosity  ratio  term  in 


Eq.  (7)  ,    ,O.H 


where ^(^   is  salt  water  viscosity  evaluated  at  tube  wall 

temperature,  is  considered  negligible  and  will  hereafter  be 

dropped  from  the  expression  of  Eq.  (7). 

Relation  (7)  is  based  upon  the  following  assumptions 

.  fully  developed  flow  in  smooth  tubes. 

.  fluid  properties  are  evaluated  at  the 
bulk  fluid  temperature. 

and  is  valid  for  the  following  condition 

J  L 

For  fully  developed  turbulent  flow  in  a  tube  as 
defined  by  Eq.  (5),  the  Dittus-Boelter  correlation  [Ref.  3] 
expressed  as 

OS       0.4 

was  used.   Nusselt  and  Prandtl  numbers,  Nuj   and  fr  ,  are 
previously  defined  by  Eq.  (9). 

Relation  (8)   is  based  upon  the  following  assump- 


tions 


.  fully  developed  flow  in  smooth  tubes. 

.  fluid  properties  are  evaluated  at  the  bulk 
fluid  temperature 


and  is  valid  for  the  following  conditions 
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.  Prandtl  numbers  ranging  from  0.6  to  100. 

.  moderate  temperature  differences  between 
the  wall  and  fluid  conditions. 

c.  Salt  Water  Fouling  Heat  Transfer  Resistance 
In  this  document,  it  will  be  assumed  that  the 

fouling  resistance  coefficient  for  tubeside  salt  water  can 
be  maintained  at  .00025  (  hr- H2  F"/37*u  )  using  one  of  the 
following  techniques: 

.  Chlorination. 

.  MAN  Brush  System. 

.  Amertap. 

.  Chemical  cleaning 
Pressure  drops  associated  with  cleaning  techniques 
will  not  be  considered  in  this  analysis.   Piping  losses  will 
be  a  function  of  tube  length,  inner  diameter,  salt  water 
velocity  and  the  absolute  roughness  of  the  tubing  design 
material  only. 

d.  Ammonia  Shellside  Heat  Transfer  Coefficient ,  h^ 
Initially,  /-j^  will  be  assumed 

kn3=  1°00  (BTU/hr-HZ-FV  (9) 

since  its  value  cannot  be  directly  calculated  during  this 
phase  of  the  analysis. 

Using  the  thermal  resistance  expressed  as 

R<   -  >.  -  .. 
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Rz    = 


R3  = 


FU  = 


da 


do   in  do  /  dj 

1 


an  initial  value  for  the  overall  heat  transfer  coefficient  may 
be  calculated. 


U0  =    (10) 


5.   NTU-effectiveness  Relations 


The  NTU-effectiveness  relationships  will  be  used  to 
determine  the  evaporator  outlet  salt  water  temperature. 
Currently,  all  salt  water  properties  have  been  based  upon 
the  initial  assumption  that 

7~BdLK   =  Tn-t   (s"W  INLET  To    e.\ii\Pj 

The  expression  for  the  number  of  transfer  units  (NTU) 
which  is  a  measure  of  the  size  of  the  heat  exchanger  is  given 
by 

NT"U  =  U0At/C„„, 

where  Cimn   ^s  defined  as  capacity  rate  of  the  single  phase 
flow  in  an  evaporative  or  condensing  two  phase  flow  regime. 
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Ctnm  -    ^5W  Cp^ 


(ID 


Evaporator  effectiveness  can  then  be  expressed  as 

C'NTU) 

£-  f-  e  (i2) 


for  two  phase  flow  regardless  of  the  flow  geometry. 
Using  the  definition  of  effectiveness 


actual  heat  transfer 

Effectiveness  =  (13) 

maximum  possible  heat  transfer 


£  -     j2_       _    ATmm    __     J~,  -  Thc  (14) 

tyfyvjix         L^.1  max  THl  —lei 

The  expression  for  /^77tfMrePresents  the  single  phase 

(salt  water)  flow  and  7c'  represents  ammonia  inlet  temperature 

to  evaporator  taken  at  state  point  3A. 

6 .   Evaporator  Salt  Water  Outlet  Temperature  and 
Bulk  Temperature 

Using  the  relationships  of  Eqs .  (12)  and  (14),  the 

following  expression  may  be  formulated  for  salt  water  outlet 

temperature 
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Concurrently,  a  revised  evaporator  average  salt 
water  temperature  can  be  expressed  as 

Tbulk    =  (Thi  +  Th,)/Z  (16) 

Using  the  revised  value  for  average  salt  water 
temperature, iterate  with  equation  (1)  until  the  revised  and 
current  values  of  bulk  temperature  satisfy  a  specified 
convergence  criterion. 

7 .   Amount  of  Heat  Absorption,  £p 

Using  the  results  of  Eq.  (16)  and  (12)  ,  the  amount  of 
heat  absorption  by  the  evaporator  may  be  expressed  as 

Q    -     Cmirl     (Tfri-    \  Ho)  (17) 

8  .   Log  Mean  Temperature  Dif  f  erence  ,  LtATQ 

The  NTU-effectiveness  method  can  be  used  to  determine 

the  mean  effective  temperature  difference  (LMfD)  across  the 

evaporator  (heat  exchanger)  . 

Using  Eq.  (17)  and  the  definition  of 

Q  --  Uo  A,  F  LtATi) 

the  log  mean  temperature  difference  across  the  evaporator 
may  be  expressed  as 
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C-NTU) 

Lnro     -   Qmnliz^ 1^-TcJ  (18) 


where  Tci    =    Isjh^   evaluated  at  state  point  3. 

F   =  correction  factor  onLMTP,  equal  to  1  for 
two  phase  flow. 
9.   Film  Temperature  for  Property  Evaluation ,  Tf 

In  order  to  evaluate  the  shellside  ammonia  heat 
transfer  coefficient,  working  fluid  properties  (i.e., 
viscosity,  specific  heat,  etc.)  must  be  evaluated  at  the  film 
temperature  to  validate  critical  heat  transfer  expressions. 

By  modifying  the  expression  in  Eq .  (10)  multiplying 
by  a  single  tube  outer  area,  a  value  for  single  tube  con- 
ductance can  be  expressed  as 

U.A    =  A 


Ri  +  Rz   r  R3  i-  fly 

Subsequently,  the  average  amount  of  heat  transferred 
per  tube  would  equate  to 


0=  U.A  (Ibulk-  h) 


where    T~3    =  Ihh*    evaluated  at  state  point  3. 

Again  using  the  resistance  analysis  in  Section  3, 
shellside  wall  temperature  may  be  derived  from 


TW2  =  Teu,K  -  q  (JL±£*-l2* 
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Knowing  shellside  wall  temperature  and  the  free-stream 
temperature,  film  temperature  can  be  derived  from  their 
arithmetic  mean. 


Z 


f 


(19) 


10.   Ammonia  Mass  Flow  Rate,  /?%>,, 

According  to  first  law  of  thermodynamics  for  steady 
state,  steady-flow  conditions  in  the  evaporator: 


EVAPORATOR  (ideal) 


Q(hot  sw) 


™HH3    h~  +  9  =  '^3  ^3 


(20) 


from  which  the  ammonia  mass  flow  rate,  YY\^5  ,  may  be  determined 
if  the  enthalpies  at  state  points  2  and  3  are  known. 

If  we  initialize  the  lower  and  upper  bounds  of  the 
analysis  in  terms  of  pressure  P,  and  P_,  respectively,  and 
initially  assume  that  a  saturated  vapor  leaves  the  evaporator, 
the  following  relations  may  be  expressed 
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113  s  n5Jp3  T~3  -  fwjp3 


(21) 


where       hi      =  represents  enthalpy  at  state  point  1  at 

the  suction  inlet  to  the  working  fluid 
circulation  pump. 
ho     =  represents  enthalpy  at  (ideal)/state  point  3 
as  a  saturated  vapor. 
llt  13    -    represent  the  respective  saturation 

temperatures . 
V-f   =  represents  the  specific  volume  at  state 
point  1. 
To  summarize,  the  upper  and  lower  pressure  bounds  of 
the  system  (  Pi  and  ?5  )  will  be  initialized  in  the  analysis 
and  treated  as  design  variables  by  the  otpimization  code. 
Temperature  at  state  point  3  is  initially  assumed  to  be  a 
saturated  vapor  (ideal  T~) ;  however,  the  working  fluid  is 
subject  to  a  shellside  pressure  drop  as  it  passes  across 
the  evaporator  with  an  outlet  quality  of  90-95%.   Properties 
at  state  point  3  (actual)  will  be  assessed  in  follow-on 
sections . 
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AMMONIA  CIRC  PUMP 


"*•——'" 


mNtt3  hi  +-  VVcf  =  wN„3  hz  (22) 

Assuming  steady  state,  steady-incompressible  flow, 
the  change  in  kinetic  and  potential  energies,  and  heat 
losses  are  negligible  for  isentropic  conditions,  and  the 
isentropic  pump  work  can  be  expressed  as 

\\'CPS     =  V1  (fz  -  Pi)  Luh&rc     Pz  -  Pss 

After  the  isentropic  pump  work  is  calculated,  the  actual 
(adiabatic)  pump  work  may  be  determined  using  pump  efficiency, 

Actual  outlet  enthalpy  at  state  point  2  may  be  deter- 
mined using  the  results  of  Eq.  (23)  with  Eq .  (22)  knowing  the 
enthalpy  at  state  point  1  from  Eq.  (21). 

Using  the  results  of  Eqs .  (21)  and  (22),  the  mass 
flow  rate  in  Eq.  (20)  may  be  calculated  as  the  average  shell- 
side  mass  flow  rate  for  the  working  fluid  (ammonia) . 
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11 .   Tube  Profile,  Flow  across  Tube  Bank,  and  Tube  Sheet 
Diameter 

Since  the  heat-exchanger  arrangements  (evaporator 
and  condenser)  involve  multiple  rows  of  tubes,  the  geometric 
arrangement  of  the  tube  profiles  is  important  in  the 
determination  of  the  heat  transfer  coefficient,  the  tube 
sheet  diameter  and  the  shell  side  pressure  drop  associated 
with  two-phase  flow  (homogeneous  model)  [Ref.  4]. 

The  following  geometric  arrangements  are  used: 


IN-LINE 


FLOW 


where  Sn 
Pp. 

A? 


=  pitch  ratio  x  outer  tube  diameter,  equal  to  Sp 
=  pitch  ratio;  the  distance  between  tube  centers 

with  respect  to  outer  tube  diameter. 
=  tube  profile  area  (centerline  to  centerline) 

per  tube . 


S/i  =  P*  da 


(24) 


Af=  Sn 


z 


(25) 
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STAGGERED 


FLOW 


5n  ■=  2  P&  d0  5  m  3C 


(26) 


SP   =  Pa  do  COS  3o 


(27) 


Therefore,  the  tube  profile  area  (centerline  to  centerline) 
per  tube  is  equal  to 


4p--  ^ 


r>.   *J3 


(28) 


The  ratio  of  minimum  flow  area  to  the  frontal  area  can  be 
expressed  as 

/iff     Si,  -  4q 


A, 


5n 


(29) 


Using  the  selected  tube  profile  geometry,  either  in- 
line or  staggered,  and  knowing  the  required  number  of  tubes 
by  equation  C2) ,  the  tube  sheet  diameter  for  heat  exchanger 
design  can  be  assessed  as  follows: 
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A 


(30) 


where  Tso    -   tube  sheet  diameter. 

To  estimate  the  shellside  ammonia  flow  velocity  the 
following  control  volume  is  introduced  (ammonia  circulation 
piping  and  the  top  portion  of  the  evaporator) . 

r  c 


'TUBE  BUNDUE.  £Bfl*iTAL  SJiREACE  -     » 

ITiiuhii  iTTTl  f 


mf 


If  the  mass  flow  rate  remains  unchanged  across  any 
boundary  (continuity) , 

Furthermore,  if  we  assume  the  evaporator  has  the 
means  to  evenly  distribute  liquid  droplets  across  the  top 
of  the  tube  bundle  (spray  nozzles  and  baffling) ,  the  follow- 
ing expressions  can  be  applied  to  estimate  the  mean  droplet 
velocity  approaching  the  bundle: 

Let  (Af.)u9   *  Af\ 

where    "^  =  percent  of  tube  frontal  area  which  is  occupied 
by  droplets. 

The  mass  flow  rates  are 
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where    A?    =  ammonia  pipe  cross -sectional  area. 

Vp  =  average  ammonia  velocity  in  the  pipe 

Therefore 

A, 


and  since 

it  follows  that  the  average  velocity  of  ammonia  through  the 
circulation  pipe  is  equivalent  to  the  average  velocity  of 
ammonia  at  the  tube  frontal  area  boundary. 

Vp^Vf  C3i) 

Thus  the  assumption  that  7(jx Ap //\f    is  equivalent  to 
the  assumption  of  constant  liquid  kinetic  energy  in  the 
transition  from  the  pipe  exit  to  the  bundle  entrance.   Con- 
sidering the  minimum  free-flow  area  for  shellside  flow  passage, 
Aff  can  be  derived  from  Eqs .  (29)  and  (30): 


Af.   =  TSj>    l~t 


Aff   =  Af    (Sn-cl.\  (32) 

>■    Oh  I 


5» 

where   Af  =  represents  the  flow  frontal  area. 
L^    =  tube  length. 
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Using  the  calculated  values  of  Eqs.  (32)  and  (20), 
the  mass  velocity  for  the  minimum  free-flow  area  can  be 
expressed 

where  })iNH  represents  the  average  ammonia  mass  flow  rate. 

12 .   Pressure  Drop  of  Two-Phase  Flow  across  a  Bank  of 
Tubes  ,  /\p 

This  portion  of  the  analysis  will  use  an  analytical 
model  for  two-phase  pressure  drops  applicable  for  a  fog  or 
spray  flow  pattern  occurring  at  high  void  fractions  --  the 
homogeneous  model  [Ref.  4]. 

The  model  asserts  that  if  the  pressure  drop  in  the 
two-phase  flow  for  a  liquid-vapor  mixture  is  relatively  small 
compared  to  the  absolute  pressure,  the  flow  is  considered 
incompressible.   Subsequently,  the  density  of  each  phase  is 
practically  constant.   During  the  process  of  phase  change, 
the  phase  and  velocity  distributions  are  changed,  and  so  is 
the  momentum  of  the  flow.   Therefore,  the  pressure  drop  of 
a  vertical  two-phase  flow  consists  of  three  components: 
friction  loss,  momentum  change,  and  elevation  pressure  drop 
arising  from  the  effects  of  the  gravitational  force  field. 

The  local  pressure  gradient  for  a  two-phase  flow 
may  be  expressed  as 

A  ?roT   =  ^PpflcnOM  *"  A  PHcMchTUV   +APelEv«tiox, 
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For  a  given  channel  length, /_^_,  the  pressure  drop 
components  can  be  represented  by 


A  P  -  £  G'  U    L 

^  '  FRICTION  '     j\    z  '  ( 


(34) 
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Hchen tu m 


9,» 
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LEvMT/OH 


vj< 


L 


c. 


and  the  total  pressure  drop ,  /^  £    ,  is  given  by  the  sum  of 
these  expressions 

where    j   =  single-phase  friction  factor  by  Jakob  expressed 
in  Eqs .  (35)  and  (36)  . 
Lc     ~   channel  flow  length,  defined  for  horizontal 
tubed  evaporators  as  Lc-  Tsj>        (tube  sheet 
diameter) . 
-De  =  equivalent  diameter  of  flow  channel,  defined 

by     As.  =  Pr  da  -  d0  . 
V"  =   mean  specific  volume  defined  by 


V 


-  %[f  +  5^jr**). 


where    X     =  quality  of  mixture  (state  point  3) . 

Vf    =  specific  volume  of  liquid  (state  point  1) 
T/G  =  specific  volume  of  vapor  (state  point  3)  . 
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The  basic  assumptions  of  the  homogeneous  model  (fog 
flow  model)  [Ref.  4]  are: 

(1)  equal  linear  velocities  of  vapor  and  liquid, 

(2)  thermodynamic  equilibrium  between  the  two 
phases,  and 

C3)  a  suitably  defined  single-phase  friction  factor 
is  applicable  to  the  two-phase  flow. 
Using  assumption  (3)  and  the  correlations  by  Jakob 
[Ref.  3] ,  a  suitable  single-phase  friction  factor  can  be 
calculated  from  previously  defined  tube  profile  relationships : 
for  staggered  tube  arrangements: 


f=      {    0.2S-+   21£5 Re  (35) 

*         1        [(Sn-iJ/do]1'*  r  M 


and  for  in-line  tube  arrangements: 

f=  {o.(m*f. o.oa sn/d.    — _   R  (36) 


[(Sn-d.)Aj0-43    ^ 


where  Reynolds  number  (max)  is  determined  from  the  shellside 
ammonia  flow  and  the  nozzling  effect  of  the  tube  geometry  as 
expressed  by 


vmax  -     Vf  ; 


where    Vc    =  the  ammonia  velocity  at  the  tube  frontal  area 
boundary  determined  by  equation  (31) . 
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Reynolds  number  for  maximum  shellside  flow  can  be 
calculated  using  the  following  expression 


Rw  >*  "T  ~  C37) 


ALf 


Eq.  (37)  and  tube  profile  data  can  then  be  used  to 
evaluate  the  single-phase  friction  factor,  required  for 
Eq.  (34).   All  other  components  of  the  total  pressure  drop 
Eq.  (33)  can  be  determined  from  previously  calculated  data. 
13 .   Pressure  Drop  Across  the  Moisture  Separator.  Apm  sap 

This  portion  of  the  analysis  will  simulate  the  use 
of  a  cyclone  separator  to  improve  te  evaporator  outlet  vapor 
quality.   The  flow  pattern  in  a  cyclone  separator  is  complex 
and  simplifying  assumptions  are  inadequate  to  allow  the 
calculation  of  the  corresponding  pressure  drop,  which  can 
vary  from  1  to  20  inlet  velocity  heads  [Ref.  5].   Therefore, 
the  worst  case  condition  will  be  applied  with  an  approxima- 
tion for  the  fluid  flow  inlet  area  to  the  separator  banks. 

By  approximating  the  inlet  area  as  a  fraction  of 
the  evaporator  frontal  area 

/Wr=  0<\  Tsp   Lt 

the  inlet  fluid  velocity  can  then  be  determined  using  the 
working  fluid  mass  flow  rate,  Eq .  (20). 


K*,*/* 


A,u.rrV 


'iNLai 


where   /£>      =  density  of  ammonia  at  state  point  3 
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Therefore,  if  the  pressure  drop  across  the  moisture 
separator  is  equal  to  20  times  the  inlet  velocity  head, 


AP  „.  =  zo/yd 


,n.scp  ■'    ?*_  (38) 


^ 


u 


14.   Enthalpy  at  State  Points  3  and  4 

Since  Eq.  (33)  represents  the  pressure  drop  across 
the  evaporator  shellside,  the  actual  pressure  at  state 
point  3  or  evaporator  outlet  may  be  determined  from 

/5(A/£WJ=  P3  -AFEVA?  (39) 

where  P3   was  previously  described  as  the  pressure  for  a 
saturated  vapor. 

Similarly  the  actual  pressure  at  state  point  4,  the 
moisture  separator  outlet,  may  be  expressed  as 

P4=  PzCNBWj-AP^sep  (40) 

Operating  under  the  dome  of  the  Temperature-Entropy 
diagram,  the  following  properties  are  defined 


Sr\3f  (new)  -  hfjpSOm)  \\4i  =  hfj  pj 

Ii3ic*&*)-  ^)^t^w)        fag*  hyP4 


(41) 


47 


The  subscript  (N£w)    representing  a  revised  property 
will  hereafter  be  dropped  from  the  expressions  in  Eq .  (41). 

Assuming  an  evaporator  outlet  quality  of  9 0 - 9 5 %  ,  and  a 
moisture  separator  outlet  quality  of  99-99.5%,  enthalpies 
at  state  points  3  and  4  may  be  determined  using  the  relation- 
ships of  Eqs .  (41) 

fl4S  k^X4C|14J-|i4-f)  (42) 

15 .   Revised  Ammonia  Mass  Flow  Rate  and  Velocity 

Till  now,  we  assumed  that  the  shellside  mass  flow 

rate  was  given  in  accordance  with  the  ideal  system  defined 

by  Eq.  (20);  however,  in  actuality  this  is  not  the  case. 

The  diagramatic  representation  that  follows  better 

illustrates  the  heat  absorption  phase  of  the  OTEC  power 

system  and  will  provide  the  basis  for  the  analysis  and 

optimization . 

Note,  as  in  the  previous  control  volume  analysis, 

the  following  conditions  are  assumed. 

.  Steady  state. 

.  Steady-incompressible  flow. 

.  Change  in  potential  and  kinetic  energies  is 
negligible . 
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Analyzing  the  moisture  separator  as  a  separate 
control  volume, 
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If  we  assume  that  there  is  no  carry-over  of  vapor  in 
the  separator  drain,  then 


and 


X3  m$   =  X4  M4 


X4  • 
X3 


(43) 


However,  for  reasons  of  flow  continuity,  the  mass 
flow  rate  through  the  separator  drain  must  be  included  in 
the  control  volume  analysis;  therefore 

fyj3  =  \V\4  ¥  mD£ 


(44) 
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Substituting  Eq .  (43)  into  (44)  and  solving  for  }^De, 
the  following  expression  can  be  derived 


X4 


rri4 


(45) 


Looking  at  the  evaporator  as  a  separate  control 


volume , 


the  energy  balance  is 


Q+  mzjnz  ¥  ^Uih^L  +  mDi  hoi  -  W3/13  rhi^iuc 

Assuming  the  change  in  enthalpy  across  the  re-flux 
pump  and  the  difference  between  the  separator  drain  outlet 
and  evaporator  inlet  are  negligible,  the  energy  balance 
becomes 


(46) 
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where   hpe  -  nf  J/P3  rfi?\   fluid  drained  from  the  separator  is 

assumed  to  be  a  saturated  liquid. 
Furthermore,  a  mass  balance  of  the  evaporator  control 
volume  can  be  expressed  as 


YY\Z   r  mm  i-moi  =  W3  f  A>1*e  C47) 

where     MA^-Wlili    . 

Solving  Eq.  (47)  for  the  mass  flow  rate  at  state 
point  3  and  substituting  into  Eq .  (46)  with  Eq .  (45)  yields 
the  following  expression 

£  f  mz  V\z  +/M.  -  j  j  nu  hoc  =  mt  ||  H3  (4g) 

In  addition,  a  mass  balance  for  steady-state,  steady-flow 
indicates  that  the  mass  flow  rates  at  state  points  2  and  4 
are  equal  and  therefore 

IV4  -  ihz  (49) 

Using  Eqs .  (48)  and  (49),  the  revised  mass  flow 
rate  at  state  point  2  may  be  determined.   Concurrently,  the 
revised  average  ammonia  velocity  acting  on  the  tube  profile 
geometry  may  be  determined  from  this  revised  mass  flow  rate. 

Using  the  revised  ammonia  velocity  acting  on  the 
tube  profile  geometry  and  iterating  from  Eq .  (31)  until 
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an  acceptable  convergence  criterion  is  achieved  provides 
the  pressure  drops  across  the  evaporator  and  moisture 
separator,  and  the  properties  at  state  points  3  and  4  for 
a  given  film  temperature.   The  result  is  more  representative 
of  the  heat  absorption  phase  in  the  OTEC  power  cycle  than  is 
the  commonly  used  ideal  analysis. 

In  addition,  solving  for  the  revised  temperature  at 
state  point  3, 

T3   =  TsatJF3  (50) 

and  iterating  through  Eq.  (18)  revises  the  film  temperature 
and  subsequent  working  fluid  properties. 

16 .   Revised  Shellside  Ammonia  Heat  Transfer  Coefficient 

In  the  search  for  acceptable  correlations  to  predict 
the  average  evaporative  heat  transfer  coefficient,  two  analyti- 
cal treatments  were  found  that  lent  themselves  to  OTEC 
power  system  conditions. 

The  first  of  these  correlations  seeks  to  predict 

thin  film  evaporation  heat  transfer  coefficient  for  horizontal 

tubes  [Ref.  6].   Owens  [Ref.  6]  uses  (1)  the  similarity 

between  evaporation  and  condensation,  (2)  the  correlation 

forms  of  local  evaporation  heat  transfer  coefficients  for 

water  on  a  vertical  tube  developed  by  Chun  and  Seban,  and 

(3)  the  dependence  of  heat  transfer  on  the  vertical  spacing 

of  the  tubes  as  was  experimentally  demonstrated  by  Liu,  to 

arrive  at  the  following  correlations  for  non-boiling  thin 

film  evaporation: 
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for  laminar  flow 


0.1  -\J5  -\/3 

h-     zz(»)  (   ««+      \   (jr.)  csi) 

U/    \3/>fzkf*J    \/4*    J 


for  turbulent  flow 

fr-  nA*s(*\  I     dd±        \    (£r*f  )  (52) 


J  ft   K93J     \    K 


where   £L   =  vertical  spacing  with  respect  to  tube  outer 

d* 

diameter. 

I    =  tube  flow  rate  per  unit  length. 

The  laminar- turbulent  transition  point  is  defined 

by  the  intersection  of  Eqs .  (51)  and  (52) 


-US' 


^  -  "*>  ( ^f ) 


f 

The  pseudo -Reynolds  number  for  horizontal  vertical  falling 
film  evaporation  is  defined  by  Ref.  7. 

Re=  l£ 

The  second  correlation  combines  boiling  and  evapora 
tion  of  liquid  films  on  horizontal  tubes,  applicable  for 
vertical  banks  of  plain  and  enhanced  tubes  [Ref.  8]. 
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The  overall  model  for  a  single  tube  is  expressed  as 

H-  Ub  -h  \u    Ld    +■  \-)c(  1-   U)  (53) 

L  ^    L  / 

where   H^  -   Rohsenow  pool  boiling  correlation  over  the 
entire  tube  length  given  by 


,   __  4^     Jk^   AT2 


Jig  1   iVjD 


3/"t 


with    Cs-f  =  function  of  the  fluid-surface  combination. 

/\  I   =  wall  temperature  minus  free  stream  saturation 

temperature . 
(Jc   =  surface  tension 
fti   =  heat  transfer  coefficient  in  the  developing 

region. 

'-  -  3  c,  r 


"'  -  i  ^ 


.4/3 


3/tf. 


f/5 


tw«-i 


47//>^    V     5/1,- 


and  he  =   fully  developed  heat  transfer  coefficient  given  for 
laminar  flow  by 

2  -1/3  -O.ZZ 
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and,  for  turbulent  flow, 

where    L  =  circumferential  length  of  heated  surface. 

c<  =  thermal  diffusivity. 

*-tf/  =  developing  length  around  tube  circumference. 

I    =  flow  rate  per  unit  axial  length  of  tube. 

To  apply  Eq.  (51)  for  a  vertical  bank  of  tubes , L 
is  expressed  as 

1-=  MtTTA/2 


The  laminar- turbulent  transition  point  is  defined  by 
the  intersection  of  Eqs .  (55)  and  (56) 

-\.D(o 

As   before,    the   pseudo-Reynolds   number   is    defined  by   Ref.    7 


Re  =    41  (57) 


After  using  Eq.  (57)  to  establish  which  flow  regime 
the  system  is  operating  in,  the  revised  heat  transfer  coeffi 
cient  for  non-boiling  thin  film  evaporation  or  nucleate 
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boiling  may  be  calculated  and  then  iterated  with  the  initial 
assumption  for  the  shellside  heat  transfer  coefficient,  Eq.  (9). 
This  will  have  a  convergence  effect  on  variables  which  are 
a  function  of  the  shellside  heat  transfer  coefficient,  moving 
them  closer  to  actual  OTEC  system  performance  characteristics. 

The  user  should  be  aware  that  the  predictions  for  the 
OTEC  power  system  using  ammonia  have  been  for  the  case  where 
no  boiling  occurs  in  the  film.   This  condition  is  dictated 
by  industrial  preference  for  plain  tube  heat  exchangers  to 
minimize  fouling  and  the  characteristic  of  ammonia  to  wet 
surfaces  well,  flooding  out  nucleation  sites.   A  number  of 
enhancement  techniques  have  been  developed  to  create 
nucleate  boiling,  including  a  variety  of  tube  configurations 
and  surface  preparations;  however,  a  preference  for  them 
has  not  materialized.   The  nucleate  boiling  development  in 
Eq.  (51)  which  would  be  indicative  of  tube  enhancement  is 
provided  for  information  only  and  will  not  be  included  in 
the  optimization  or  summary  of  conclusions. 

Having  described  the  methods  used  to  predict  the 
shellside  heat  transfer  coefficient,  we  can  complete  this 
chapter  of  the  OTEC  power  system  analysis  by  constructing  the 
heat  exchanger  cost  analysis. 
17 .   Evaporator  Cost  Analysis 

At  the  request  of  TRW,  Wyatt  Industries,  a  large 
exchanger  fabricator,  prepared  cost  estimates  for  three 
different  sizes  of  vertically  configured  evaporators  and 
condensers,  based  upon  initial  design  specifications  prepared 
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by  TRW.   Based  upon  these  estimates,  TRW  developed  sets  of 
equations  that  represent  the  costs  of  various  heat  exchanger 
component  parts  for  shell  diameters  ranging  from  10-35  ft 
and  35-50  ft  [Ref.  9] . 

The  following  are  the  TRW  evaporator  cost  ($) 
equations  as  a  function  of  outer  tube  diameter  (inch) , 
total  number  of  tubes  and  tube-sheet  diameter  (ft)  for 
tube-sheet  diameters  of  10-35  ft. 

.  Drilling  time/tube  sheet  thickness 

td=  o.u  (da-  o.s)  (58) 

.  Thickness  of  the  tube  sheet 

0,63 

±rs  =  0>Sh    l*o  (59) 

.  Tube  sheet  labor  cost 
dTSL^1SGCf^(N±/U30)(i:el/o.uX±rs/^)  (60) 

.  Tube  sheet  material  cost 

Cm*  I**-** LTs?  {61) 

.  Tube  installation  cost 

o.l 
Cn=  34  Ni  do  (62) 

.  Heat  exchanger  shell  cost 
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.  Ammonia  distribution  plate  and  battles  cost 
CWB  =  l38bS.7*(Kt/lt,30Xtd/0.u)CTS0  /iS)  *  (64) 

.  Bustle,  flanges  channels  and  flow  plates  cost 

CBFCF   =  308SSO  (T^/1S)Z  (65) 

.  Tube  material  cost 

CTH,  (e\  Lt  +  Ez)  N±  A  (66) 

where    ti  =  curve  fit  of  tube  cost  per  foot. 
£2   =  tube  machining  cost  if  required 

.  Heat  exchanger  head  costs 

Cm=  53240  (T,o/lsJ  (67) 

.  Water  inlet,  nozzles  and  supports  cost 

Cwis=22COI0,k(Ts,/16)  (68) 

.  Tube  welding  costs  (Titanium  tubes) 
for     hit    -  3hCCC 


1.03  ,       i  O*  7 

Ct*j  =  14'73  W*  (do /is) 


(69) 


for  AL  >360PO 


m  A4-'  (d./\.f) 

The  sum  of  cost  Eqs .  (60)  through  (69)  would  equal 
the  cost  to  fabricate  one  OTEC  evaporator  with  a  tube  sheet 
diameter  of  10-35  feet  (all  the  preceding  component  costs 
have  been  adjusted  for  current  pricing  at  a  10%  annual  rate 
of  inflation) . 
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If  our  analysis  is  based  on  a  30-year  life-cycle 
criterion,  no  adjustments  are  necessary  to  any  component 
cost  equation  if  titanium  tubing  is  used  due  to  its  anti- 
corrosive  qualities;  however,  using  aluminum  tubing  (i.e., 
Al-5052),  the  expense  of  retubing  must  be  considered  to  meet 
the  criterion  of  a  10-year  life  cycle  for  aluminum  tubing. 
This  implies  Eq.  (61)  and  (66)  must  be  modified  to  reflect 
the  costs  of  retubing  at  the  10  and  20-year  point  in  the  cycle. 

.  Aluminum  tube  installation  cost 

Cwx  a  £rr  [  1  «■  O  t  i/V  (  f  ♦  </°]  (70) 

where    i        =  projected  inflationary  rate  (input  by  customer) 
.  Aluminum  tube  material  cost 


^^[ifO^Ao^n 


(71) 


For  tube  sheet  diameters  of  35-50  ft  the  following  cost 

relationships  apply  [Ref .  9] : 

.  Equations  for  drilling  time/tube  sheet 

thickness  (58) ,  thickness  of  tube  sheet  (59) , 
and  tube  material  costs  remain  unchanged. 

.  Tube  sheet  labor  and  material  cost  (titanium) 

0,-Wi       0>(o8 
tT5L--S$.^^  Tsp         td  (72) 


2,014 
C'r5Af  =  Z^'SU    r*P         id  (73) 
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Tube    sheet   labor   and  material    cost    (aluminum) 
CT5L=  73.181  Mb  Tso         ±d  (74) 


1.4,1 

CT5H^  Z&.sTsd         tT5  (75) 


Tube    installation   costs 

0.1 
C-l  ~  3k.5~42  AJt  do  (76) 


Heat  exchanger  shell  cost 

C,»;--^^(L*^^  (77) 


Ammonia  distribution  plate  and  baffle  costs 

1.82  0,873 

CWB-  \SB.tXlT5t       +  12.411  Ut       hd  (78) 

Bustle,    flanges,    channels,    and   flow  plate   costs 

_     2.1Z 

CbFCF  =4iz.<m  fsp  (79) 

Heat   exchanger   head  cost 

C-HXH    =    mS'5]    Ts»  C80) 


Water   inlet,   nozzles   and   support   cost 

Cvtms  =  144S-ZW  TidU  C81) 
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Tube   welding   costs    (titanium   tubes) 
for       j\Jt    4:  Z&COG 

Crw=  Hl3Nt'°\clo/i.s)  C82) 


for      A/+  >  5&0OO 

1,03/        .        \0>1 

Cr^  =  0.811-7  N±      (do/i.s) 


As  indicated  previously,  the  cost  to  fabricate  one 
OTEC  evaporator  with  a  tube  sheet  diameter  35  to  50  ft  is 
equal  to  the  sum  of  component  costs  Eqs .  (72)  through  (83) 
(all  the  preceding  component  costs  have  been  adjusted  for 
current  pricing  at  a  10%  annual  rate  of  inflation) . 

For  an  analysis  based  on  a  30-year  system  life-cycle 
criterion,  the  additional  costs  for  aluminum  retubing  must 
be  considered  and  Eqs.  (70)  and  (71)  apply. 
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IV.   PARASITIC  LOSSES 

A.   INTRODUCTION 

This  chapter  describes  in  detail  the  programming  analysis 
for  parasitic  losses  which  include:   (1)  pumping  and  pipe 
requirements  for  both  cold  and  hot  salt  water  systems, 
(2)  pumping  and  pipe  requirements  for  the  working  fluid 
(ammonia)  circulation  and  re-flux  systems,  and  (3)  turbine 
generator  losses  due  to  inefficiencies.   Hotel  requirements 
have  not  been  incorporated  into  the  analysis,  but  could  be 
included  for  the  final  design  analysis. 

Pumping  power  requirements  will  be  determined  through 
the  use  of  the  general  energy  equation  between  the  inlet  and 
outlet  of  the  system  control  volume  [Ref.  3]. 

/L  z  z 

dt+Vi     +   a  Zi    =     Vo    +.  a  H„  +  \{Js  +  (Losses) 

0 

To  determine  the  pumping  power  Ws   the  following  effects 
will  be  evaluated: 

1.  Density  head. 

2.  Friction  losses. 
.  Intake  piping. 

.  Heat  exchanger  tubing. 

.  Exit  piping  (if  employed) . 

3.  Thermodynamic  pressure  head. 

4.  Elevation  head. 
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5.   Minor  losses. 

.  Intake  piping  inlet  configuration  (contraction)  . 

.  Intake  piping  screen  (obstruction) . 

.  Flow  through  valves,  elbows,  etc. 

.  Outlet  piping  (expansion) . 

.  Inlet  to  heat  exchanger  tubing  (contraction) . 

.  Outlet  from  heat  exchanger  tubing  (expansion) . 

.  Outlet  of  exit  piping  (if  employed) . 

In  the  above  pump  head  evaluations,  the  following  inputs 

are  specified: 

.  Pipe  lengths  (hot,  cold,  ammonia  circulation  and 
re-flux  piping) . 

.  Inner  pipe  diameters  (initialized  and  treated  as  a 
design  variable  by  the  optimization  code) . 

.  Absolute  roughness  corresponding  to  piping/tubing 
material  (designer  specified) . 

.  Fluid  velocity  (initialized  and  treated  as  a 
design  variable  by  the  optimization  code) . 

.  Pump  mechanical  and  electrical  efficiencies. 

As  an  overview  of  the  parasitic  pump  loss  analysis,  the 

following  major  steps  in  the  algorithm  are  listed  in  order 

of  their  execution: 

.  Hot  pipe  salt  water  pump. 

. .  Inlet  piping  friction  losses  (86)  . 

. .  Minor  piping  losses  due  to  inlet  screen  (87) 
and  plenum  design  to  evaporator  core  (88). 

..  Evaporator  core  minor  losses  (89,  90)  and 
tubeside  friction  losses. 

. .  Total  pressure  losses  (92)  and  pumping 
head  (93) . 

..  Pumping  power  requirements  (95). 
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. .  Pump  cost  analysis  (96)  . 

Cold  pipe  salt  water  pump. 

. .  Initialize  cold  pipe  inner  diameter  and  SW 
velocity  (design  variables) . 

. .  Minor  losses  due  to  inlet  ducting  (97)  and 
plenum  design  to  condenser  core  (98) . 

. .  Inlet  piping  friction  losses  (99) . 

..  Condenser  core  minor  losses  (100,  101)  and 
tubeside  friction  losses  (103)  . 

. .  Density  head  (104) . 

..  Total  pressure  losses  (105). 

. .  Pumping  power  requirement  (107) . 

..  Pump  cost  analysis  (108). 

Ammonia  circulation  pump. 

. .  Piping  friction  (109)  and  minor  losses  due  to 
valving/elbows  (110) . 

. .  Pressure  drop  across  evaporator  shellside  (112) 

. .  Thermodynamic  head  (113) . 

. .  Elevation  head  (114) . 

.  .  Total  pressure  losses  (115) . 

. .  Pumping  power  requirement  (116)  . 

. .  Pump  cost  analysis  (118) . 

Ammonia  re -flux  pump. 

. .  Piping  friction  (119)  and  minor  losses  due  to 
valving/elbows  (120)  . 

. .  Thermodynamic  head  due  to  pressure  drop  of 
saturated  liquid  ammonia  across  evaporator 
shellside  (122). 

..  Elevation  head  (123). 

..  Total  pressure  losses  (124). 
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. .  Pumping  power  requirements  (126) . 
..  Pump  cost  analysis  (127). 
.  Parasitic  pump  losses. 
In  the  following  section,  the  basic  steps  summarized  above 
will  be  described  in  detail. 

B.   ANALYSIS  OF  PARASITIC  LOSSES 

1 .   Hot  Pipe  Salt  Water  Pump,  PHP 

The  pressure  losses  due  to  piping  friction  and 
associated  minor  losses  will  be  determined  using  the 
Darcy-Weisbach  correlation  [Ref.  10]. 


n  z 

KlV 


1=1  /     zfr 


(83) 


where   fcl   describes  the  resistance  coefficient. 
V   =  fluid  velocity. 

where    -f      =  friction  factor. 

—   =  equivalent  length  in  pipe  diameters. 

P 

In  order  to  determine  the  friction  factor,  the  pipe 

flow  Reynolds  number  must  be  calculated. 


(84) 


p       /^S-tO   Vsa)    oil 
^~   ~ 


yU 
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where  A\^,Ms^~   properties  of  wait  water  at  the  hot  pipe 

inlet  temperature  (assumed  constant  through- 
out the  pipe)  . 
V<  i  d'L    =  salt  water  velocity  and  inner  pipe  diameter 
(initialized  and  treated  as  design  variables 
by  the  optimization  code) ;  velocity  assumed 
constant  over  pipe  length. 
Pipe  flow  Reynolds  number  greater  than  2300  will  be 
considered  turbulent. 

for  laminar  flow 


(85) 


for  turbulent  flow 


i    - 


[ln(£/3.7rfi  ^H/fct/^] 


(86) 


where  €   =  absolute  roughness  corresponding  to  piping  material 
selected. 
Eq.  (86)  yields  a  friction  factor  within  one  percent 
of  the  Colebrook  equation  and  is  valid  for  the  following 
conditions  [Ref.  9], 

D 

a 
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Considering  the  resistance  coefficient  for  pipe  minor 

losses 

.  Assume  the  inlet  duct  is  the  same  size  as  the 
pipe  inner  diameter,  but  it  is  screened 

K=f-5"  (87) 

.  Assume  piping  enters  evaporator  through  an 
area  which  is  abruptly  changed  [Ref.  11] 

z 
K=  [l-  (Ji/Tsof]  (88) 

where   Tsd     ~   evaporator  tube  sheet  diameter  (assume  tube 

sheet  diameter  is  twice  as  large  as  the  inner 
pipe  diameter) . 

Summing  the  results  of  Eqs .  (84),  (87),  and  (88)  to 
determine  the  total  resistance  coefficient,  the  pressure 
losses  due  to  piping  can  then  be  determined  using  Eq.  (83). 

If  a  variety  of  valves  or  fittings  are  to  be  included 
with  Eq.  (84),  Ref.  11  provides  a  representative  listing  of 
equivalent  length-to-pipe-diameter  values. 

To  analyze  the  pressure  drop  across  the  evaporator 

tubeside,  we  again  use  the  Darcy-Weisbach  correlation,  but 

for  different  design  assumptions. 

Assume  inlets  to  evaporator  tubing  are  well 
rounded  [Ref.  11] 


K=  cr 


(89) 
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Assume  outlets  of  evaporator  tubing  expand 
to  an  infinite  reservoir  [Ref.  10] 

<=  I*0  (90) 


Using  the  Reynolds  number  in  the  previous  chapter, 
Eq.  (5),  the  corresponding  friction  factor  Eq.  (85)  or  (86), 
and  resistance  coefficient  can  be  determined 


where    L^.,  d't      =  evaporator  tube  length  and  inner  tube 

diameter  and  are  initialized  and  treated  as 
design  variables  by  the  optimization  code. 
Summing  the  results  of  the  resistance  coefficient 

in  Eqs .  (89),  (90)  and  (91),  the  pressure  losses  due  to  the 

evaporator  design  may  be  determined  using  the  Darcy-Weisbach 

correlation  Eq.  (83). 

The  results  of  the  piping  losses  and  core  design 

losses  are  equivalent  to  the  hot  pipe  salt  water  pumping 

system  requirements 

A??amp    =  APpiP^^EM    ^A^EVAP  D£SI**J  (9  2) 

converting   to   pumping   head 

H-  \     APpun,?  (93) 
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Pumping  power  in  terms  of  horsepower  can  be  determined 
using  the  following  expression 


p    _    msw  (  6  H 


3 


(94) 


where   "Tfa  =  pump  mechanical  efficiency  (designer  input)  . 

/f]5w  =  salt  water  mass  flow  rate  determined  in 
previous  chapter,  Eq.  (2). 

To  equate  parasitic  pump  losses  to  power  input, 
Eq.  (94)  is  converted  to  the  motor  load  requirement  in  terms 
of  megawatts  electrical. 

tnP^j-  LOLx    Cch^E-iisiok)   FACToK  (95) 

where   "7Lu*  =  pump  motor  efficiency  (designer  input)  . 

Because  of  the  high  salt  water  flow  rates  and  rela- 
tively low  pumping  heads,  good  engineering  design  would 
dictate  the  use  of  axial  flow  (propeller)  type  pumps. 

Using  the  algorithm  developed  by  TRW  [Ref.  9]  from 
data  provided  by  Johnston  Pump  Co.,  and  Process  Equipment 
Co.  (distributors  of  Ingersoll  Rank  and  Johnston  Pumps), 
the  cost  of  salt  water  pumps  can  be  expressed  as 


C 


(b/iooo)o,7o  +Scl  1.2U  to  (96) 
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where 


D=  ii£i"v. 

4 


z. 


where   d[    \/2^     =  inner  hot  pipe  diameter,  salt  water  velocity 

(initialized  for  analysis  and  treated  as 
design  variables  by  the  optimization  code) . 
The  above  algorithm  is  valid  for  the  following  con- 
ditions 

.  vertical,  wet  pit,  propeller  type  pumps  with 
cast  iron  steel  columns  with  protective  epoxy 
coating,  stainless  steel  shaft  and  bronze 
impeller. 

.  pump  size  from  155,000  through  750,000  GPM 
with  total  dynamic  heads  of  8  through  12  feet. 

Eq.  (96)  has  been  adjusted  for  current  pricing  at  a  101 

annual  rate  of  inflation. 

2.   Cold  Pipe  Salt  Water  Pump, ftp 

Using  Reynolds  number 

to.  /?*»}   Vr*)  ch 

U  ~  yUsvJ 

where  /Qw }M^    -   properties  of  salt  water  at  the  cold  pipe 

inlet  temperature  (assumed  constant  through- 
out the  pipe) . 
Vsfij    d'L    ~   salt  water  velocity  and  inner  pipe  diameter 
(initialized  and  treated  as  design  variable 
by  the  optimization  code) ,  velocity  assumed 
constant  over  pipe  length. 
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Pipe  flow  characteristics  and  friction  factor  can  be  identi- 
fied.  A  pumping  analysis  will  be  developed  for  the  cold 
pipe  pump  using  the  Darcy-Weisbach  correlation,  similar  to 
the  development  in  the  preceding  section. 

Considering  the  resistance  coefficient  for  minor  pipe 
losses 

.  Assume  the  inlet  duct  is  well  rounded  [Ref.  11]. 


K  ~„Qa*  (97) 


INi.C-1 


.  Assume  piping  enters  condenser  through  an 
area  which  is  abruptly  changed  [Ref.  10] . 


K  =  f  i  «■  (  di  iTso) 


Z-* 


pLc-rfam   L.  J 

where  \50   =  condenser  tube  sheet  diameter  (assume  tube  sheet 

diameter  is  twice  as  large  as  the  inner  pipe 
diameter) . 


Assume  one  ninety-degree  elbow  is  required 
in  system  [Ref.  11] . 


(98) 


L  -  30 


Summing  the  results  of  Eqs .  (84),  (97),  and  (98), 
the  total  resistance  coefficient  can  be  expressed  as 

*•*(•%+ iy****1- **<*">        (99) 
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where   Lp   =  length  of  cold  pipe. 

cAl     =  inner  diameter  of  cold  pipe. 

Pressure  losses  due  to  piping  can  then  be  determined 
using  the  Darcy-Weisbach ,  Eq.  (83). 

In  analyzing  the  pressure  drop  across  the  condenser 

tubeside,  the  Darcy-Weisbach  correlation  is  used  again,  but 

for  different  design  assumptions. 

.  Assume  inlets  to  evaporator  tubing  are  well 
rounded. 

K=°'S  (100) 

.  Assume  outlet  of  condenser  tubing  expands  to 
an  infinite  reservoir. 

K=/.e?  (ioi) 

Defining  Reynolds  number  for  condenser  tubeside 
flow,  while  assuming 

I  Bulk   -    I  CoLoii^Ler)  (102) 


t\ed  -    -f-- ■ 


where  /^^    /Us\tJ  ~   properties  evaluated  at  condenser  tubeside 

bulk  temperature  (initially  assumed  equal  to 
cold  pipe  inlet  temperature) . 
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VsW,<^=  average  salt  water  velocity  through 

tubing,  inner  condenser  tube  diameter  (both 
are  initialized  and  treated  as  design 
variables  by  the  optimization  code) . 

The  corresponding  friction  factor,  Eq.  (85)  or  (86),  and 

resistance  coefficient  can  be  determined 

Kaxe.  =  f  -*  (103) 

dl 

where  Lt,d[,  the  condenser  tube  length  and  inner  tube  diameter 
are  initialized  and  treated  as  design  variables  by  the  optimi- 
zation code. 

Summing  the  results  of  the  resistance  coefficient 
in  Eqs .  (100),  (101),  and  (103),  the  pressure  losses  due  to 
the  condenser  design  may  be  determined  using  the  Darcy-Weisbach 
correlation,  Eq.  (83). 

A  complete  analysis  of  cold  pipe  losses  must  also 
include  the  effect  of  density  head  and  a  corresponding  increase 
in  pumping  power  requirements. 

For  most  engineering  problems  involving  the  flow  of 
liquids  through  a  pipe,  where  the  temperature  change  in  the 
pipe  is  small,  the  density  of  the  fluid  is  considered  to  be 
a  constant  and  the  fluid  is  termed  "incompressible."  However, 
the  flow  problem  in  OTEC  cold  pipe  systems  is  unique.   We  can 
continue  to  assume  that  there  is  negligible  change  in  the  fluid 
temperature,  virtually  unaffected  by  the  ocean  thermal 
gradients,  because  of  the  system's  characteristic  high  mass 
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flow  rates.  However,  the  height  of  the  water  column  (1500 
to  3000  feet)  inside  the  pipe  requires  the  effect  of  fluid 
compressibility  to  be  taken  into  consideration. 

The  effect  of  an  increase  in  density  with  depth 


can  be  expressed  by  the  integral 

MP 

f>3 


i 

2 


with  a   density  head  defined  as 


y\a    -    ife      —  l 


-r 


\i? 


Integrating  the  pressure-density  variation,  the 
density  head  reduces  to  [Ref.  12] 

where   ICW  ■  mean  compressibility  of  wait  water,  f (salinity, 
temperature  and  pressure) . 

fi>      =  reference  density  at  which  K^,  is  evaluated. 

Considering  pressure  at  any  depth  obtained  from  the 
integral , 


p-  -  4  (  />(z)d-2. 


the  density  head  can  be  rewritten  as  follows 


f- 


*J  I  Z 


2Note  that  2   is  measured  as  positive  upward  so  that  ocean 
depth  values  (ze_/  2t )    are  negative  and  (7l&.-2l)   is  a  positive 
quantity. 


Rigorous  procedures  for  calculating  the  density 
profile  which  is  a  function  of  temperature,  salinity  and 
pressure  may  be  found  in  Ref.  13;  however,  they  will  not 
be  discussed  in  this  document. 

For  the  purposes  of  simplification,  the  following 
solution  technique  was  developed: 

(1)  If  the  liquid  in  the  pipe  is  taken  to  have  a 
constant  density  with  respect  to  pressure,  the  compressibility 
approaches  zero;  the  density  head  can  then  be  expressed  as 

(2)  Converting  the  geometric  term  for  elevation  to 
an  equivalent  integral  expression 

The  reference  density  is  taken  to  be  the  inlet  value  so  that 


and  the  density  head  can  be  rewritten  as  follows 


2- 1 

(3)  Assuming  a  linear  distribution  of  density  with 
depth,  due  to  temperature  variations,  as  illustrated  below 
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2e         Zi 
the  following  linear  expression  for  density  with  respect  to 
depth  may  be  formulated,  where  "Z^r^      for  convenience. 


yOi  ~/0=(/>\  -f>c)(  I  "  */Zi) 


C4)  Applying  the  equation  developed  in  section  3  to 
the  density  head  integral  above  and  integrating  over  the 
range  of  values  for  sea  water  depth  (z) ,  the  following 
equation  is  derived  as  a  linear  approximation  to  the  density 
variation  of  sea  water  with  respect  to  depth 


where    Z?L  ,/?<?  =  curve  fit  evaluations  of  density  for 

specified  depths  of  sea  water.   Data 
extracted  from  Ref.  14. 
The  results  of  the  piping  losses,  core  design  losses, 

and  density  head  are  equivalent  to  the  cold  pipe  salt  water 

pumping  system  requirements 


Zji  pu/np  =  ZjrJ 


pipe.  SLjsreM 


AP 


Cduo  DE-Stfti 


+■ 


AP 


fretisiry 


(105) 
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Using  Eq.  (93),  Eq.  (105)  can  be  converted  to  a 
pumping  head.   Similarly,  pumping  power  in  terms  of  horse- 
power can  be  determined  using  Eq.  (44) . 


where 


z 


(106) 


and 


^yvJ   =  density  of  salt  water  evaluated  for  a 

constant  inlet  temperature. 
Vs\u  ai    =  cold  pipe  salt  water  velocity,  and  inner 

diameter  (initialized  and  treated  as  design 
variables  by  the  optimization  code) .   Note 
salt  water  velocity  through  cold  pipe  is 
considered  to  be  constant. 
Pumping  power  can  then  be  expressed  in  terms  of 
megawatts  electrical 

hlP^vo;  =  LIZ   X  C0rfV£K*IO,\J  FACTOR  (10  7) 

where  ~tfH=   pump  motor  efficiency  (designer  input) . 

Using  the  same  arguments  for  the  selection  of  an 
axial  flow  (impeller  type)  pump,  as  used  for  the  hot  pipe 
salt  water  pump,  the  pump  cost  algorithm  developed  by  TRW 
can  be  applied  to  the  cold  pipe  salt  water  pump  assuming 
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the  required  conditions  are  validated. 

(  P/lOOo)o.1£  r^olj.ZUlD  (108) 


C 


Equation  (108)  has  been  adjusted  for  current  pricing  at  a 
101  annual  rate  of  inflation. 

3.   Ammonia  Circulation  Pump^  Fcirc 

The  function  of  the  ammonia  circulation  pump  is  to 
circulate  and  lift  saturated  liquid  ammonia  from  the  condenser 
hot  well  at  state  point  1  and  increase  its  pressure  to  exceed 
the  operating  conditions  in  the  evaporator  at  state  point  2. 

In  order  to  evaluate  these  characteristics,  the 
following  pumping  elements  will  be  included  in  the  analysis: 
.  Piping  losses  (friction  and  minor) . 
.  Heat  exchanger  shellside  pressure  drop. 
.  Thermodynamic  pressure  head. 
.  Elevation  head. 
As  in  the  preceding  analysis,  Reynolds  number  is  used 
to  determine  pipe  flow  characteristics 


Re.,  =  A 


yd, 


where  /?)M-    -   saturated  liquid  properties  of  ammonia  for 

the  temperature  at  state  point  1  (assume  any 
temperature  increase  from  pump  work  is  negligible) 
di    =  inner  pipe  diameter  (initialized  and  treated  as 
a  design  variable  by  the  optimization  code) . 
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V  =  ammonia  flow  velocity  determined  from  the  pre- 
ceding chapter,  Eq.  (50). 
The  ammonia  pipe  friction  factor  can  then  be  deter- 
mined from  Eqs .  (85)  or  (86),  and  the  piping  friction 
resistance  coefficient  can  be  expressed  as 

d, 

where  L  =  ammonia  circulation  pipe  length  (designer  input) . 

Considering  the  resistance  coefficient  for  minor  pipe 
losses,  assume  there  are  four  ninety-degree  elbows  in  the 
system 

K=4i  (110) 

where  —  =*  equivalent  length  in  pipe  diameters  for  a  standard 
elbow  [Ref.  11]. 
Summing  the  results  of  Eqs.  (109)  and  (110),  piping 
losses  (friction  and  minor)  can  be  determined  using  the  Darcy- 
Weisbach  equation  (83). 

APoo-=A *-  *4-A]-^  (in) 


The  heat  exchanger  shellside  pressure  drop  is  also 
included  in  the  pumping  head  requirement  because  it  serves  as 
a  resistance  to  flow. 
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Pressure  drop  across  the  evaporator  shellside  was 
determined  using  the  two-phase  flow  model  (homogeneous) 
expressed  by  Eq.  (33) 

^-'evap  ~  ^fpRlCT/oM  *&  MOHZtiruH'*' &**£L£VArtOM  (112) 

Since  the  pump  is  required  to  lift  the  working  fluid 
to  a  higher  elevation  and  increase  its  operating  pressure, 
the  following  elements  must  be  included  in  the  analysis: 
.  Thermodynamic  head 

where      AR  =  R>  -  Pi  (113) 

*-3  THERMO 

represents  the  difference  in  thermodynamic 
operating  pressure  between  state  point  2  and 
state  point  1. 
.  Elevation  head 

where  APEl£VATl0li  =  ?Z  "  21  (114) 

'£\    =   datum. 

^Zz    -  elevation  of  the  evaporator  inlet  above 
datum  (taken  to  be  equal  to  evaporator 
tube  sheet  diameter  plus  25) . 
represents  the  lift  head  required  to  move  the 
working  fluid  to  a  higher  elevation. 
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The  results  of  piping  losses  (111) ,  evaporator 
pressure  drop  (112)  ,  the  thermodynamic  head  (113)  and 
elevation  head  (114)  are  equivalent  to  the  ammonia  circulation 
pump  system  requireements . 

APpun)p~  ^M ?IP£   +  ^i£VAP*  A*  THERMO  ^  A  t  ELEVATION  (115) 

Using  Eq.  (93)  with  ammonia  properties,  Eq.  (115) 
can  be  converted  to  pumping  head  and  finally  expressed  as 
pumping  power  (horsepower) . 

a--  %*  (-#) 

■ 

where     /^nh3  =  mass  flow  rate  of  ammonia  determined  by 

Eq.  (20)  of  the  previous  chapter. 
"Tjp   =  pump  mechanical  efficiency  (designer  input) . 
Pumping  power  can  then  be  expressed  in  terms  of  mega- 
watts electrical 


fc\KHMw)'      —X  CjD*IUe*SlOM    FACTOR  (117) 


where  *^|  =  pump  motor  efficiency  (designer  input)  . 

Because  of  high  pumping  head  and  moderate  flow  rates, 
good  engineering  design  would  dictate  the  use  of  a  single 
suction  centrifugal  flow  type  pump. 
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Using  the  algorithm  developed  by  Westinghouse  Electric 
Co.  [Ref.  15]  from  data  provided  by  Bingham  Pump  Division, 
Portland,  Oregon,  the  cost  of  the  ammonia  circulation  pump 
can  be  expressed  as 


CpUr»F=     ("V-?  It)        1.21*10  (118) 

1  30100     J 


where  fr\kla_    =  mass  flow  rate  of  ammonia  (^bm/hr) 

Vr      =  specific  volume  of  saturated  liquid  ammonia 
at  state  point  1  (ft  /\bm) 
Eq.  (118)  has  been  adjusted  for  current  pricing  at  a  10% 
annual  rate  of  inflation. 

4.   Ammonia  Re -flux  Pump,  Pre-flux 


The  function  of  the  re-flux  pump  is  to  recycle 
ammonia  droplets  which  are  not  evaporated  in  the  heat  absorp- 
tion process.   Saturated  liquid  at  approximately  the  heat 
exchanger's  operating  pressure  is  lifted  from  the  evaporator 
drain  to  the  ammonia  feed  inlet,  for  redistribution  as  droplets 
across  the  evaporator  tube  bundle.   (Drainage  mass  flow  rate 
is  assumed  to  be  equal  to  30%  of  the  evaporator  inlet  feed 
mass  flow  rate.) 

In  order  to  evaluate  these  characteristics,  the 
following  pump  elements  will  be  analyzed: 

.  Piping  losses  (friction  and  minor) . 
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.  Thermodynamic  pressure  head. 
.  Elevation  head. 
As  in  the  preceding  analysis,  Reynolds  number  is  used 
to  determine  pipe  flow  characteristics 


■J 


R*,  --  ^ 


where  /^iM   =   saturated  liquid  properties  of  ammonia  for  the 
average  pressure  across  the  evaporator. 
di      -    inner  pipe  diameter  (initialized  and  treated 

as  a  design  variable  by  the  optimization  code), 
v   =  ammonia  flow  velocity  determined  from  the 
evaporator  drainage  mass  flow  rate  assumed 
equal  to  301  of  the  evaporator  inlet  feed  mass 
flow  rate  (assume  velocity  constant  throughout 
the  pipe) . 
The  re-flux  pipe  friction  factor  can  be  determined 
from  Eqs .  (85)  or  (86),  and  the  piping  resistance  coefficient 
can  be  expressed  as 

K~  f—  (119) 

di 

where  L  ~   ammonia  re-flux  pipe  length  (designer  input) . 

Once  again,  considering  the  resistance  coefficient 
for  minor  pipe  losses  assume  there  are  four  ninety-degree 
elbows  in  the  system 
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K.4± 


(120) 


where  L.   =  equivalent  length  in  pipe  diameters  from  a  standard 
D 

elbow. 

Summing  the  results  of  Eqs .  (119)  and  (120),  piping 

losses  (friction  and  minor)  can  be  determined  using  the 

Darcy-Weisbach,  equation  (83) 


AP 


piPc 


L  (X,  ^  J 


(121) 


In  order  to  determine  the  thermodynamic  pressure 
head,  the  pressure  drop  across  the  evaporator  for  the 
saturated  ammonia  liquid  must  be  analyzed.   Since  the 
saturated  vapor  and  liquid  are  in  thermodynamic  equilibrium, 
the  results  of  Eq.  (112)  apply.   Therefore 

APu^P*-Pz. 

Therefore,  the  thermodynamic  pressure  head  is  equal 
to  the  pressure  drop  across  the  evaporator  for  the  saturated 
ammonia  liquid. 


A  P 


TtteRMO 


=  APU 


9 


(122) 


Finally,  the  elevation  head  is  equal  to  the  elevation 
of  the  evaporator  feed  inlet  with  respect  to  datum,  the 
drain  outlet. 
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Therefore , 

4PBW-  £z-£\  (123) 


ELSV 


where    Z\    =  datum,  drain  outlet. 

^£  =  elevation  of  the  evaporator  inlet  above  datum 
(taken  to  be  equal  to  the  evaporator  tube 
sheet  diameter  plus  10) . 
The  results  of  piping  losses  (121)  ,  the  thermodynamic 
pressure  head  (122)  ,  and  elevation  head  (1233  are  equivalent 
to  the  ammonia  re- flux  pump  system  requirements. 


AP       =/\P     +-AP  +  4P  t124) 

^'pamp         '-1    pipe.  BLBVATiOM        ^    THE&MO 


As  before,  using  Eq.  (93),  Eq.  (124)  can  be  converted 
to  a  pump  head  and  finally  expressed  in  terms  of  pumping 
power  (horsepower) . 

fiU*.»!a  (*h)  (125) 

where   /T)a   =  drainage  mass  flow  rate. 

y\p     =  pump  mechanical  efficiency  (designer  input) . 
Pumping  power  can  be  expressed  in  terms  of  megawatts 
electrical 

fi*S-FU^MwF  ^^*  CON\Z£d5IOisJ  FACTOK  U26) 
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where  7lH   =  pump  motor  efficiency  (designer  input) . 

Using  the  same  arguments  for  the  selection  of  a 
centrifugal  pump,  the  pump  cost  algorithm  developed  by 
Westinghouse  can  also  be  applied  to  the  ammonia  re-flux 
pump. 

c^r-($^)  w*">  cm) 

s 

where   1'Y}^    =  mass  flow  rate  of  evaporator  drainage  ammonia 

Ur     =  specific  volume  evaluated  at  the  average 

evaporator  pressure  (ft / \bm) 
Eq.  [127)  has  been  adjusted  for  current  pricing  at 
a  10%  annual  rate  of  inflation. 
5 .   Parasitic  Pump  Losses 

Parasitic  pump  losses  is  the  summation  of  electrical 
auxiliary  pumping  requirements  (hotel  and  maintenance  loads 
not  included)  determined  by  Eqs .  (95),  (107),  and  (126). 

?loss   =    Fh?   +  Per    <"  Po«c  <-  Prs-fiux  <128) 
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V.   TURBINE  AND  ELECTRICAL  POWER 

A.   INTRODUCTION 

The  turbine  generator  is  one  of  the  critical  elements 
of  the  OTEC  power  system.   Its  energy  conversion  efficiency 
and  efficiency  of  design  have  a  major  effect  on  the  overall 
system  performance.   To  illustrate  this  point,  Ref.  16 
reported  that  a  three-point  change  in  turbine  efficiency  from 
85  to  88%  results  in  a  3.6°o  increase  in  gross  power,  and  a  5% 
increase  in  net  power  developed. 

This  chapter  will  describe  the  analysis  to  evaluate  the 
expansion  turbine  thermodynamic  properties  and  generator 
output.   The  use  of  these  properties  will  determine  the 
internal  turbine  efficiency  and  outlet  quality  subject  to 
design  and  thermodynamic  constraints.   The  relationship 
between  the  condenser  operating  pressure  (design  variable) 
and  the  turbine  outlet  quality  will  be  used  to  initialize 
the  heat  rejection  characteristics  of  the  condenser. 

General  literature  on  turbomachinery  designed  for  OTEC 
closed  cycle  systems  indicates  that  a  turbine  having  the 
following  characteristics 

.  Double  flow,  axial  inflow, 

.  Four  stages  of  expansion, 

.  Operating  at  1800  RPM, 
provides  the  optimum  aerodynamic  design  [Ref.  16].   However, 
it  is  not  the  intent  of  this  thesis  to  analyze  the  geometry 
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and  performance  parameters  of  the  turbine.   Turbine  geometry 
such  as 

.  Specific  speed  and  specific  diameter, 

.  Wheel  diameter, 

.  Rotational  speed, 

.  Blade  height, 

.  Blade  stresses, 
should  be  treated  as  a  separate  systems  problem  using 
optimization  to  improve  state-of-the-art  design. 

Parasitic  losses  due  to  the  following  generator  turbine 
inefficiencies  will  be  evaluated  in  this  section. 

.  Generator  mechanical  and  electrical. 

.  Turbine  mechanical. 

As  an  overview  of  the  turbine-generator  analysis,  the 

following  major  steps  of  the  algorithm  are  listed  in  order 

of  their  execution: 

.  Gross  electrical  output  with  no  parasitic  losses 
(129). 

.  Enthalpy  at  state  point  5  (130) . 

.  Turbine  outlet  quality  (131)  . 

.  Entropy  at  state  point  from  a  specified  outlet 
quality  (132)  . 

.  Quality  and  enthalpy  at  state  point  5s  (133,  134). 

.  Internal  (adiabatic)  turbine  efficiency  (135) . 

.  Turbine  cost  analysis  (137)  . 

.  Generator  cost  analysis  (138)  . 

In  the  following  section,  the  basic  steps  summarized  above 

will  be  described  in  detail. 
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B.   ANALYSIS  OF  THE  TURBINE  AND  ELECTRICAL  POWER  REQUIREMENTS 

1 .  Gross  Electrical  Output  and  Inefficiency  Losses 

If  the  net  electrical  output  required  is  indicated 
by  (in  terms  of  megawatts) ,  the  gross  electrical  load  at 
the  turbine  shaft  can  be  expressed  as 

E^iTyf    *  Ploss  (129) 

where    Ploss       =  parasitic  pump  losses  determined  by  Eq .  (128) 
?1tm   =  turbine  mechanical  efficiency  (designer 

input) . 
y[a£N    ~   generator  mechanical  and  electrical 
efficiency  (designer  input) . 
The  loss  of  electrical  output  due  to  generator- 
turbine  inefficiencies  is  equal  to 

E      =E/ — ^ ) 

2 .  Turbine  Efficiency 

The  power  developed  across  the  turbine  is 

E^  =  m  (/?*•- In) 

where      M       =  mass  flow  rate  of  ammonia  given  by  Eq.  (48). 

r)4  =  enthalpy  at  state  point  4,  Eq.  (42). 
From  this,  the  enthalpy  at  state  point  5  can  be  calculated. 
If  we  initialize  the  operating  pressure  of  the  con- 
denser in  terms  of  r5 ,  the  following  relations  may  be  expressed 
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^5  =  ^yp>  h$+   =  /if  J  P^  (130) 

Therefore,  it  follows  that  the  turbine  outlet  quality, 
XS",  can  be  determined  from 


hs  =  hi"f  +  X5"(/l5J-  /tffj 


(131) 


Having  established  the  moisture  separator  outlet 
pressure  and  temperature,  Eqs .  (40)  and  (41),  the  entropy 
at  state  point  4  can  be  determined  for  a  known  separator 
outlet  quality  (designer  input)  using  the  following 
relations 

54  =  S*f  v-  X4  (  54J  -  54f )  (132) 

For  isentropic  turbine  work, 

54  -  Sss  (133) 

the  quality  at  state  point  5s  may  be  determined  using  the 
following  relations 

5o9  -  $<■]  Its'  5i'f  =  S+-  J  TS* 


5  r  ^3^5- 

Ss*  =  Srf  +  X55-(5^  -  5^1         (134) 

Having  determined  the  quality  at  state  point  5s,  the 
enthalpy  can  now  be  determined. 
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hfs  =  Jns-f  +  Xi~5  (Jnsy  -  hof) 


(135) 


Using  the  results  of  Eqs .  (41),  (130),  and  (132), 
the  internal  turbine  efficiency  (adiabatic)  can  be  deter- 
mined, expressed  by 


v\r'  ~i       i 


r"    j       ,  (136) 

v    P,4  -  YiSS 

To  ensure  a  realistic  selection  of  internal  efficiency, 
the  following  constraints  are  attached  to  the  optimization 
code 

#    Ins-  <  fi5~5 

'  X5'5<XS' 

3.   Turbine  Cost  Analysis 

The  ammonia  turbine  cost  is  based  on  an  algorithm 
developed  by  Westinghouse  to  estimate  manufacturing  costs 
[Ref.  15] . 

<W  M**(KS15-+L/{3l.i>oollt)Ff.  (137) 

where    E^  =  gross  electrical  output  in  KW. 
Aj'^  =  2  (for  a  double  -flow  turbine)  . 
fr     =  flow  price  factor  (1.0  for  single-flow,  1.447 
for  double-flow) . 
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The  above  algorithm  is  valid  for  the  following 
conditions : 

.  Double  flow,  axial  inflow. 
.  Multi-stage. 
.  Operating  at  1800  RPM. 
The  generator  cost  will  be  based  on  an  algorithm 
developed  by  TRW  from  data  provided  by  selected  manufacturers, 

C^£a/-    (p'023Ef+0*3j1.2\Xlb  (138) 

and  is  valid  for  the  following  conditions 

.  1800  RPM  rotor  speed. 

.  power  factor  0.8. 
Eqs .  (137)  and  (138)  have  been  adjusted  for 
current  pricing  at  a  10%  annual  rate  of  inflation. 
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VI .   CONDENSER 

A.   INTRODUCTION 

As  indicated  in  the  introduction  to  Chapter  III,  several 
heat  exchanger  concepts  have  been  proposed  for  the  closed- 
cycle  OTEC  system,  with  variations  in  their  design. 

The  analysis  to  be  presented  for  the  condensing  heat 

exchanger  will  be  based  upon  the  following  design 

characteristics : 

.  Single-pass  shell  and  tube  heat  exchanger. 

.  Horizontal/vertical  orientation  of  tubes  with  an 
equilateral  triangle  or  square  tube  profile. 

.  Smooth  plain-tube  configuration  (no  enhancements) . 

.  Tube  material  (titanium  or  aluminum  based  on  a 
30-year  life-cycle  criterion). 

.  Biofouling  control  based  upon  an  achievable 
fouling  factor. 

.  Heat  exchanger  centerline  located  on  sea  surface. 

As  an  overview  of  the  condenser  analysis,  the  following 

major  steps  in  the  algorithm  are  listed  in  order  of  their 

execution : 

.  Initialization  of  design  variables  (DV) . 

. .  Tube  length. 

..  SW  velocity  through  condenser  tubes. 

..  Outer  tube  diameter. 

..  Tube  profile  pitch  ratio. 

.  Amount  of  heat  rejection  (139) . 

.  Tubeside  bulk  temperature  (142)  . 
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.  Total  number  of  tubes  (143) . 

.  Log  mean  temperature  difference  (144)  . 

.  Conductance  (146) . 

.  Number  of  transfer  units  (145)  . 

.  Heat  exchanger  effectiveness  (147) . 

.  Initially  assume  a  value  for  ammonia  heat 
transfer  coefficient  (151) . 

.  Single  tube  conductance  (148)  . 

.  Average  heat  rejection  per  tube  (152). 

.  Film  temperature  (153)  . 

.  Revised  ammonia  heat  transfer  coefficient  (154,  etc.) 
iterate  with  (151) . 

.  Tube  profile,  flow  parameters  across  the  tube 
bank  (158  ,  etc.)  . 

.  Tube  sheet  diameter  (163) . 

.  Condenser  shellside  pressure  drop  for  two-phase 
flow  (166)  . 

.  Revised  properties  at  state  point  1  (171,  172); 
iterate  with  (21)  . 

.  Overall  heat  transfer  coefficient  (173)  . 

.  Total  heat  transfer  surface  area  (174)  . 

.  Revised  condenser  tube  length  (175)  . 

.  Heat  exchanger  cost  analysis. 

In  the  following  section,  the  basic  steps  summarized 

above  will  be  described  in  detail. 

B.   ANALYSIS  OF  THE  CONDENSER 

« 

1 .   Amount  of  Heat  Rejection, Q ' 

Using  the  calculated  value  for  enthalpy  at  state 
point  5,  equation  (131)  from  the  previous  chapter,  the  ideal 
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values  at  state  point  1,  Eq.  (21),  and  the  steady-state  mass 
flow  rate  of  ammonia,  Eq.  (48),  the  amount  of  heat  rejected 
by  the  condenser  can  be  expressed  as 

9=  'nw,  (hff-Ju)  (139) 

2 .   Tubeside  Bulk  Temperature 

As  in  condenser  tubeside  Reynolds  number,  salt 
water  properties  will  be  evaluated  at  bulk  temperature, 
initially  assumed  equal  to  the  cold  pipe  inlet  temperature. 

Using  this  premise,  the  condenser  salt  water  capacity 
rate  can  be  evaluated 

d^/^mcpcp^  d40) 

where  Cp       =    specific  heat  of  salt  water  initially  evaluated 
at  the  cold  pipe  inlet  temperature. 

i 

jfyi<in   =  mass  flow  rate  of  salt  water  through  the  cold 
pipe  previously  evaluated  by  Eq.  (107). 

Using  the  results  of  Eqs .  (139)  and  (140),  and  the 
known  cold  pipe  inlet  temperature,  the  condenser  salt  water 
outlet  temperature  may  be  evaluated  from  the  basic  expression 

Cy-  Cm,n  (Tc0-TcJ  (i4i) 

where    7c*  Tci    ~   condenser  salt  water  outlet  and  inlet 

temperatures,  respectively. 
Having  determined  the  condenser  salt  water  outlet 
temperature,  the  revised  bulk  temperature  can  be  expressed  as 
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Tx=    Tc0*  '*i  C142) 


Using  the  revised  condenser  bulk  temperature  and 
iterating  with  Eq.  (102)  corrects  the  operating  temperature 
for  salt  water  properties  which  are  essential  to  the 
analysis . 

3.   Total  Number  of  Condenser  Tubes,  Kit 

Since  the  mass  flow  rate  of  salt  water  through  the 
cold  pipe  is  equivalent  to  the  mass  flow  rate  through  the 
condenser,  according  to  the  law  of  continuity, 

it  follows  that  the  number  of  condenser  tubes  for  a  specified 
tube  diameter,  can  be  evaluated  using  the  following  expression 


4- 

Elys™  Tdi  V^Hj  f143) 


4 

where  /^sw    ~   average  salt  water  density  evaluated  at 
bulk  temperature. 
di     ■  inner  tube  diameter  (initialized  and  treated 

as  a  design  variable  by  the  optimization  code) . 
Vt  =  average  salt  water  velocity  through  the  con- 
denser (initialized  and  treated  as  a  design 
variable  by  the  optimization  code) . 
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4 .   Log  Mean  Temperature  Difference, LMTD 

Using  the  result  of  Eq.  (141),  the  known  pipe  salt 
water  inlet  temperature,  and  the  inlet  temperature  of 
ammonia  evaluated  at  state  point  5,  the  LMJ~D  of  the  condenser 
may  be  expressed  as 


5 .   NTU-Effectiveness  Relations 

The  number  of  transfer  units  which  is  a  measure  of 
the  condenser  size  can  be  determined  from  the  basic 
expression 

WU  =  U«  A.  (145) 

where  the  conductance  (Ual\A  of  the  heat  exchanger  is  a 
function  of  the  heat  absorbed  and  the  LMTD . 

Q=(li,A„)LMrD  C146) 

The  condenser  effectiveness  can  then  be  expressed  as 

(-MTU) 

<£-)-e  (147) 

for  a  two-phase  flow,  regardless  of  the  flow  geometry. 
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6  .   Single-Tube  Conductance,  LL/\, 

Using  the  resistance  analysis  derived  in  Chapter  III, 
Section  4  for  an  initialized  tube  length 

L=L; 

the  heat  exchanger  conductance  for  a  single  tube  can  be 
expressed  as 

1 


LI.  A.- 


_i_    +    -L'R         |  IncUdl     ,  J-p    |    !       (148) 


where   r]jrw  =  tubeside  heat  transfer  coefficient. 

Rfsvi    =  salt  water  fouling  heat  transfer  resistance. 

K     -   thermal  conductivity  of  the  tube  material. 
Atf.Aj  =  total  outer  and  inner  tube  surface  areas 

(including  fin  and  bare  tube) ;  tube  length 
is  initialized  and  treated  as  a  design 
variable  by  the  optimization  code) . 
KfA/ri  =  ammonia  fouling  heat  transfer  resistance 
~Ho,y\    ~   outer  and  inner  total  fin  efficiency 
a.   Tubeside  Reynolds  Number 

Since  the  salt  water  heat  transfer  correlation 
is  dependent  on  tubeside  flow,  Reynolds  number  must  be 

evaluated         _     si      \/      J- 

rfej=  ' ; 

4  /As* 

where   /^iw  Msw~   salt  water  density  and  viscosity  are 

evaluated  for  the  fluid's  bulk  temperature. 
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di^Vs*!   =  inner  diameter  and  average  salt  water  tube 

velocity. 

Reynolds  numbers  greater  than  2300  will  be 

indicative  of  turbulent  flow  [Ref .  3] . 

b.   Salt  Water  Heat  Transfer  Coefficient,  n^ 

Once  again  the  empirical  relationship  proposed 

by  Sieder  and  Tate  [Ref.  3]  will  be  used  for  laminar  heat 

transfer  in  tubes  and  as  defined  by 

1/3    J/3    0.14 

W*.*«(M.)(£J(&) 

Nusselt  and  Prandtl  numbers  are  defined  as 


hla,= 


_     hat/  d\. 


d-    k^r  (149) 


pr=  cftw  Ms*  [1S0} 


where  dynamic  viscosity,  specific  heat,  and  thermal  conduc- 
tivity of  salt  water  are  evaluated  at  the  salt  water  bulk 
temperature . 

The  effect  of  the  viscosity  ratio  in  the 
Sieder-Tate  equation  is  considered  negligible,  and  will 
hereafter  be  dropped  from  the  expression.   The  assumptions 
and  validity  condition  associated  with  the  Sieder-Tate 
equation  were  stated  in  Chapter  III,  Section  4,  and  will  not 
be  repeated  here. 
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For  fully  developed  turbulent  flow,  again  the 
Dittus=Boelter  correlation  [Ref.  3]  was  used. 

Muja  D.CZ3  R<tf'8pr°'5 

Nusselt  and  Prandtl  numbers  are  previously 
defined  by  Eqs .  (149)  and  (150).   Assumptions  and  conditions 
for  validity  were  stated  in  Chapter  III,  Section  4. 

c.  Salt  Water  Fouling  Heat  Transfer  Resistance 
As  indicated  previously,  it  will  be 

assumed  that  the  foulding  resistance  for  tubeside  salt 
water  can  be  maintained  at  . 00025 (Hr -H 'F/BTu) 

d.  Ammonia  Shellside  Heat  Transfer  Coefficient ,  h^H 
Initially,  fo    will  be  assumed 

huu    =   iCOO  (BTU/hh-H*F)  (151) 

since  its  value  cannot  be  directly  calculated  during  this 
phase  of  the  analysis. 

Using  the  following  single-tube  thermal 
resistance  i 


In  do/ dx 


Rs=      zr'KL 


Rr-  < 


\o  h,H3Td0  L 


100 


an  initial  value  for  single  tube  conductance  (outer  tube 
surface)  may  be  calculated 


\~Jio  r\o   ~ 


7.   Film  Temperature  for  Property  Evaluation,  7f 

In  order  to  evaluate  the  shellside  ammonia  heat 
transfer  coefficient,  working  fluid  properties  must  be 
evaluated  at  the  film  temperature. 

This  can  be  accomplished  by  using  the  results  of  the 
single  tube  conductance,  the  tube  side  bulk  temperature  and 
the  working  fluid  saturation  temperature,  expressed  in  the 
following  equation  for  single  tube  heat  transfer  rate  (average) 

Q=  UA0(Ts- -TBliLK)  (152) 

Again  using  the  resistance  analysis  as  in  Chapter  III, 
the  shellside  wall  temperature  may  be  expressed  as 

Tvi'TeaLK+GiRi+Rz  +  Rs) 

Knowing  the  shellside  wall  temperature  and  the  free- 
stream  temperature,  the  film  temperature  can  be  derived  from 
their  arithmetic  mean 

Tf  =  T*j2  r  Tsr  (153) 

Z 

For  purposes  of  this  calculation,  saturated  tempera- 
ture conditions  at  state  point  5  are  taken  to  represent 
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free-stream  conditions,  when  in  fact  the  two-phase  process 
will  experience  a  pressure  drop  and  a  corresponding  drop  in 
temperature . 

8 .   Revised  Shellside  Ammonia  Heat  Transfer  Coefficient, 

This  analysis  will  include  correlations  for  both 
horizontal  and  vertical  heat  exchangers. 

In  the  horizontal-tubed  condenser,  Nusselt's 
correlation  was  used  as  a  predictor  [Refs.  7  and  17], 
for  laminar  flow 

\f/3 

hxMff  %L?±JL±r)  C154) 

where       10       =  estimate  of  ammonia  mass  flow  rate  across 

each  tube. 
Kf>/V/-/f  =  properties  evaluated  at  film  temperature. 
L   =  tube  length  (initialized  and  treated  as 

a  design  variable  by  the  optimization  code) . 
This  correlation  is  probably  conservative,  since  it 
does  not  consider  turbulence  due  to  high  vapor  velocity 
or  splashing  of  condensate  [Ref .  7] . 

For  turbulent  flow,  Nusselt's  correlation  is  increased 
by  10%  as  recommended  by  Jakob  [Ref.  17] 

1/3 

(155) 


£=  1045"  I  <£jf£jl) 
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The  laminar- turbulent  transition  point  is  defined  by 
a  Reynolds  number  of  2100,  where  the  pseudo-Reynolds  number 
for  film-type  condensation  on  horizontal  tubes  is  defined 
as  [Ref.  7] 


/?«  = 


.  21 


Mh 


where  I  -  mass  flow  rate  of  condensate  per  tube  over  its 
length. 

In  the  vertical  tubed  condenser,  both  Nusselt's  and 
Kirkbride's  correlations  were  used  as  predictors  [Ref.  7]. 

For  laminar  flow,  Nusselt's  correlation  is  increased 
by  a  factor  of  1.28  as  recommended  by  McAdams  [Ref.  7]: 


h  =  U8 


\A1 


AU 


Kfft    * 


)  Mir) 


C156) 


where  J""1-  mass  flow  rate  of  condensate  per  tube  over  its 
diameter. 
For  turbulent  flow,  Kirkbride's  correlation  is  applied 

'1/3  ,C4 


\\s  o.ooii 


Kf/>*3 


Ms 


C157) 


The  laminar-turbulent  transition  point  is  defined  by 
a  Reynolds  number  of  1800,  where  the  pseudo -Reynolds  number 
for  film-type  condensation  on  vertical  tubes  is  defined  as 
[Ref.  7] 
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After  using  the  pseudo -Reynolds  number  to  establish 

the  flow  in  which  regime  the  system  is  operating,  the  revised 

heat  transfer  coefficient  for  film-type  condensation  may  be 

calculated  and  then  iterated  with  the  initial  assumption 

for  the  shellside  heat  transfer  coefficient,  Eq .  (151). 

Once  again  this  will  have  a  convergence  effect  on  variables 

in  which  the  shellside  heat  transfer  coefficient  is  a 

function,  moving  closer  to  actual  OTEC  system  operating 

point  characteristics. 

9.   Tube  Profile,  Flow  across  Tube  Bank,  and  Tube 
Sheet  Diameter        ------ 

Since  the  condenser  tube  bundle  involves  multiple 
rows  of  tubes,  the  geometry  of  the  tube  profile  arrangement 
is  important  to  determine  the  shellside  heat  transfer  coeffi- 
cient, the  tube  sheet  diameter  and  the  shellside  pressure 
drop  associated  with  the  "homogenous"  two-phase  flow  model 
[Ref.  4]. 

Using  the  same  arrangements  shown  in  Chapter  III, 
Section  2, 


IN-LINE 


FLOW 


5n  =  Pa  <*i 


AP  =  5n 


(158) 
(159) 
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where 


5« 


pitch  ratio  x  outer  tube  diameter, 
pitch  ratio  (initialized  and  treated  as  a 
design  variable  for  the  optimization  code) 
tube  profile  area  per  tube 

k-sp  — 


STAGGERED 


where 


FLOW 


5n-  2  P«  do  sin  50 


5p  =  Pr  ci0  cos  30 


the  ratio  of  minimum  flow  area  to  the  frontal  area  can  be 
expressed  as 


(160) 
(161) 
(162) 


Af 


s« 


(163) 


Using  the  selected  tube  profile  geometry  and  know- 
ing the  number  of  condenser  tubes  by  Eq.  (143),  the  tube  sheet 
diameter  for  the  condenser  design  can  be  evaluated  from  the 
following  expression 
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<vtAP^inv 


(164) 


where   Tsp  =  Tube  sheet  diameter. 

To  analyze  the  shellside  ammonia  flow  velocity,  the 
following  control  volume  is  introduced  (turbine  generator 
discharge  and  top  portion  of  the  condenser) . 

.._Js.._, 


'TUBE  BUNDLE  FRStLTAlSUHEME.-,- J 

[i  n  n  n.rTnTn  r 


Since  the  mass  flow  rate  remains  unchanged  across 
any  boundary 

Ms-  =  Mf 

Furthermore,  if  we  assume  the  condenser  has  the 
capability  to  evenly  distribute  vapor  across  the  tube  bundle 
(distribution  baffles),  the  following  development  applies 
to  the  vapor  coverage: 


Let 


where    77 


foL,  = A^ 


i 


percent  of  tube  frontal  area  which  is  covered 
by  vapor. 


ms  =  /?-  As  Vs- 
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where   /\^-    =  condenser  inlet  cross -sectional  area 
\/V  =  turbine  discharge  ammonia  velocity. 


Therefore 


If  yi^/ZrAs/A-f/^.  •>  it  follows  that  the  turbine  discharge 
velocity  is  equal  to  the  average  velocity  of  ammonia  at  the 
tube  frontal  area  boundary.   A  determination  of  the  distri- 
bution fraction  77  requires  a  detailed  knowledge  of  the 
design  of  the  turbine/condenser  interface.   In  the  absence 
of  this  information  it  is  assumed  that 

Vf  =  Vr 

A  similar  argument  could  be  presented  for  a  vertical 
tubed  condenser  where  turbine  discharge  is  admitted  to  a 
distribution  ring  that  bands  the  condenser  tube  bank. 
Exhaust  vapor  would  travel  radially  through  the  tube  bundle 
and  then  collect  at  the  bottom  after  vertical  film- 


condensation. 


Vside 


I i_. 

Again,  in  the  absence  of  a  detailed  design,  it  is  assumed 

that 


VW  -  l/V 
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Considering  the  minimum  free-flow  area  for  a  horizontal 
tubed  condenser,  A\^-f    can  be  derived  using  Eq.  (163)  and  the 
projected  frontal  area. 


where    Af    =  the  flow  frontal  area. 

Lt  =  tube  length. 
For  vertical  condensers 

fi^r  -  TTTso  *  Frontal  lz^th  of  vapor  isjle.'  FLd\aj 
Using  the  previously  calculated  value  of  the  ammonia 
flow  rate  and  Eq.  (165),  mass  velocity  for  the  minimum  free 
flow  area  can  be  expressed  as 


C=   —  (166) 


10 .   Pressure  Drop  of  Two-Phase  Flow  across  a  Bank  of 
Tubes,/\p 

The  pressure  drop  in  the  two-phase  flow  condensing 
heat  exchanger  will  be  determined  using  the  homogeneous 
model  introduced  in  Chapter  III.   The  model  will  consist  of 
three  components  --  friction  loss,  momentum  change,  and 
elevation  pressure  drop  arising  from  the  effects  of  gravity. 

The  local  pressure  drop  for  a  two-phase  flow  may  be 
expressed  as 
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AP      =  AP  MP  +-AP 


(167) 


For  a  given  channel  length,  Lc  >    the  pressure  drop 
components  can  be  expressed  by 


AR 


•f  g'v 


MOMENTUM  Qc 


u 


AP 


3 


ELEVATION  y  qc 


L 


(168) 


(169) 


(170) 


where   -f 


a 


XT  = 


single  phase  friction  factor  by  Jakob  expressed 

in  Eq.  (35)  or  (36) . 

mass  flow  velocity  determined  from  Eq.  (166). 

channel  flow  length,  defined  for  horizontal 

tubed  condensers  as  Lc-  >sd      (tube  sheet 

diameter)  and  for  vertical  tubed  condensers 

as  Lc-L-t     (tube  length). 

equivalent  diameter  of  flow  channel,  defined  by 

j)e  =  Y\  do  —  do 

mean  specific  volume  defined  by 


V«  -v. 


'■f 


O^(^-h); 
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where    X   =  quality  of  mixture  (state  point  5) . 

Vj.    =  specific  volume  of  liquid  (state  point  1) 
'Va    =    specific  volume  of  vapor  (state  point  5) . 
All  components  of  the  pressure  drop  model  Eqs .  (168,  169, 
and  170)  can  be  determined  using  the  preceding  information. 
11 .   Revised  Properties  at  State  Point  1 

Since  Eq.  (167)  represents  the  pressure  drop  across 
the  condenser  shellside,  the  actual  pressure  at  state  point  1 
or  condenser  outlet  may  be  determined  from 

where  pf  is  previously  described  as  the  condenser  operating 
pressure  for  the  ideal  cycle. 

Operating  on  the  saturated  liquid  line  on  the 
Temperature -Entropy  diagram,  the  following  properties 
are  defined: 

hl(W£w;=  r\i]pUH£w)      TimwJ-TsATjpifazw)  (172) 

The  subscript  ^AJcwJ  representing  a  revised  property 
will  hereafter  be  dropped  from  the  expression  in  Eq.  (172) . 

Until  now,  we  assumed  the  condenser  outlet  tempera- 
ture and  pressure  were  designed  to  operate  as  an  ideal  system, 
without  a  pressure  drop.   Therefore,  using  the  revised 
temperature  at  state  point  1  and  iterating  over  the  range 
from  Eq .  (21)  until  an  acceptable  convergence  criterion  is 
achieved,  all  the  preceding  variables  as  function  of  11 
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will  be  reevaluated  to  complete  the  closed-loop  cycle  of  the 
simulated  OTEC  power  system. 

12 .   Overall  Heat  Transfer  Coefficient ,  LL 

The  quantity  "U"  represents  a  measure  of  the  total 
thermal  resistances  in  the  flow  path.   Therefore,  using  the 
tube  conductance  expressed  in  Eq .  (148)  which  is  divided  by 
the  outer  heat  transfer  surface  area  of  a  single  tube,  the 
overall  heat  transfer  coefficient  for  the  condenser  can  be 
determined. 

The  thermal  resistances  are  now  expressed  as 


R<  = 


fo« 


do 


P_3  =    &o  \r\  do/  dj 


R,~-       j 


*[-  hm 


and  the  overall  heat  transfer  coefficient  for  the  condenser 
may  be  calculated  using 


R\  f  R2  t-  Rz  +Rs  t1") 
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13 .  Total  Condenser  Heat  Transfer  Surface  Area,  Ai 
Having  determined  the  corrected  number  of  condenser 

transfer  units  (145)  ,  salt  water  capacity  rate  (140)  and 
overall  heat  transfer  rate  (173)  ,  the  total  condenser  heat 
transfer  area  can  be  calculated  from  the  NTll  expression 

mti  i  -   LLAt 

N'  U        7  (174) 

14 .  Revised  Condenser  Tube  Length 

Using  the  total  heat  transfer  surface  area  calcu- 
lated from  Eq.  (174)  and  the  total  number  of  condenser 
tubes  (143)  ,  the  revised  condenser  tube  length  can  be  deter- 
mined from  the  basic  expression 

A^M^Kdo  Lt  [Revise)  (175) 

At  this  time,  it  is  necessary  ot  iterate  the  condenser 
design  until  the  two  values  (initial  and  revised)  of  the 
tube  length  converge.   This  iteration  may  be  accomplished  by 
the  COPES  routine  if  the  following  constraint  is  defined 


Lpiff   -  Lt  (Reviseo)-  Lt(r^'/r^0 


Minimization  of  this  difference  will  cause  continual 
adjustment  of  the  required  tube  length,  already  treated  as 
a  design  variable  by  the  optimization  code. 
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15.   Condenser  Heat  Exchanger  Cost  Analysis 

As  indicated  in  Chapter  III,  TRW  developed  sets  of 

equations  to  represent  the  costs  of  various  heat  exchanger 

component  parts  for  shall  diameters  ranging  from  10-35  feet 

and  35-50  feet  [Ref .  9] . 

The  following  are  the  TRW  component  cost  equations 

for  the  condensing  heat  exchanger.   Prior  equation  reference 

numbers  will  be  substituted  where  equalities  exist  with  the 

evaporative  heat  exchanger  component  cost  expressions. 

for  tube  sheet  diameter  10-35  feet 

.  Drilling  time/tube  sheet  thickness.  (58) 

.  Thickness  of  the  tube  sheet.  (59) 

.  Tube  sheet  labor  cost.  (60) 

.  Tube  sheet  material  cost.  (61) 

.  Tube  installation  cost.  (62) 

.  Heat  exchanger  drill  cost.  (63) 

.  Ammonia  distribution  plate  and  baffles 
cost . 

.  Bustle,  flanges,  channels  and  flow 
plate  cost. 

zo 

LBFCF=  UtS't&t*  Tsp  (177) 

.    Tube   material    cost. 

CrH=  (a  Lt+£t)N±  dc/itr  (178) 
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where  £l  =  curve  fit  of  tube  material  cost 
per  foot. 
(Ll>  =   tube  machining  cost  if  required. 
.  Heat  exchanger  header  cost.  (67) 

.  Water  inlet,  nozzles  and  support  cost. 

Cmns-1010(,.41S'Tsd  (179) 

.  Tube  welding  costs  (Titanium  tubes) .         (69) 
The  sum  of  the  preceding  costs  would  equal  the  cost 
to  fabricate  one  OTEC  condenser  with  a  tube  sheet  diameter 
of  10-35  feet  (all  the  preceding  component  costs  have  been 
adjusted  for  current  pricing  at  a  101  annual  rate  of 
inflation) . 

If  our  analysis  is  based  on  a  30-year  life-cycle 
criterion,  no  adjustments  are  necessary  to  any  component  cost 
equation  if  titanium  tubing  is  selected.   However,  using 
Al  5052-0,  the  expense  of  retubing  must  be  considered  to  meet 
the  30-year  life-cycle  criterion,  as  in  the  cast  of  the 
evaporation.   For  convenience,  and  possible  subsequent 
modification,  these  considerations  are  repeated  here. 

Based  upon  the  utility  of  Al  5052-0,  two  complete 
condenser  retubings  will  be  required  to  meet  the  basic  30-year 
criterion.   This  implies  Eqs .  (62)  and  (178)  must  be  modified 
to  reflect  the  costs  of  retubing  at  the  10  and  20  year  point 
in  the  cycle. 
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.  Aluminum  tube  installation  cost. 

_.  \0       .  .  20-, 

^rr  =  ^rr[H0^  <"0W  J         (180> 

where  L  =  projected  annual  inflationary 
rate  (input  by  customer) . 
.  Aluminum  tube  material  cost. 

Cm~-  £TM[\tl\  +  i)    r(\tL)    J  (181) 

for  tube  sheet  diameter  35-50  feet. 

.  Drilling  time/tube  sheet  thickness  (58) 

.  Thickness  of  the  tube  sheet.  (59) 

.  Tube  sheet  labor  and  material  costs 

(titanium).  (72,  73) 

.  Tube  sheet  labor  and  material  costs 

(aluminum).  (74,  75) 

.  Tube  installation  cost.  (76) 

.  Tube  material  cost.  (178) 

.  Heat  exchanger  shell  cost.  (77) 

.  Ammonia  distribution  plate  and 
baffles  cost. 


C„«  >  i.&zs  Nt      U  (182) 


Bustle,  flanges,  channels  and 
flow  plate. 

CSKF  =  382.824  TV2'""  C183) 


Heat  exchange  head  cost 

1.43 


(184) 
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.  Water  inlet,  nozzles,  and  supporters 
cost. 

i.OS'C 
£iviMs  =  14S3.&  T50  (185) 

.  Tube  welding  cost  (titanium  tubes) .  (82) 

As  indicated  previously,  the  cost  to  fabricate  one 
OTEC  condenser  with  a  tube  sheet  diameter  of  35  to  50  feet 
is  equal  to  the  sum  of  component  costs  (note,  all  the 
preceding  component  costs  have  been  adjusted  for  current 
pricing  at  a  10%  annual  inflation  rate)  . 

For  an  analysis  based  on  a  30-year  life-cycle 
criterion,  the  additional  costs  for  replacing  aluminum 
tubing  must  be  considered  and  Eqs .  (180)  and  (181)  apply. 
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VII.   NUMERICAL  OPTIMIZATION 

A.   INTRODUCTION 

Nearly  all  design  processes  attempt  the  minimization  or 
maximization  of  some  parameter  or  design  objective.   For 
the  design  to  be  acceptable,  it  must  satisfy  a  set  of  con- 
straints which  impose  limits  or  bounds  on  design  parameters. 

For  the  stated  problem  a  computer  program  can  be  written 
to  perform  the  basic  analysis  of  the  proposed  design.   If 
any  parameters  fall  outside  the  prescribed  bounds,  the  design 
engineer  changes  the  parameters  and  re-runs  the  program.   In 
effect,  the  computer  code  provides  the  analysis  with  the 
engineering  making  the  actual  design  decisions. 

A  logical  extension  to  the  computer-aided  approach  is  a 
fully  automated  design,  where  the  computer  also  makes  the 
actual  design  decisions  and  performs  trade-off  studies.   The 
COPES  program  provides  this  automated  design  and  trade-off 
capability  by  the  use  of  the  optimization  program  COPES/CONMIN 
[Ref.  18].   COPES  is  an  acronym  for  Control  Program  for 
Engineering  Synthesis,  and  CONMIN  is  an  acronym  for  CONs trained 
function  MINimization.   Subsequently,  a  FORTRAN  analysis 
program  simulating  a  closed-cycle  OTEC  power  system  can  be 
coupled  to  the  COPES  program  for  automated  design,  using 
some  basic  programming  guidelines  [Ref.  18] . 
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B.   COPES/CONMIN 

There  are  many  numerical  optimization  schemes  available 
to  the  engineer.   Methods  employed  by  these  schemes  fall 
into  four  basic  categories:   random  search,  sequential 
unconstrained  minimization,  optimality  criteria,  and  direct 
constrained  optimization.   The  optimization  program, 
selected  for  automated  design  analysis  of  the  simulated 
OTEC  power  system,  is  based  upon  direct  constrained 
optimization. 

Before  any  discussion  of  the  optimization  technique, 
basic  definitions  are  summarized  for  convenient  reference 
[Ref.  19] : 

.  Design  variables  -  those  parameters  which  the 

optimization  program  is  per- 
mitted to  change  in  order  to 
improve  the  program. 
.  Objective  function  -  the  parameter  which  is  to  be 

minimized  or  maximized  during 
the  optimization  process. 
.  Inequality  constraint  -  one-sided  conditions  which 

must  be  satisfied  for  an 
acceptable  design. 
.  Equality  constraint  -  condition  which  must  be 

equaled  for  the  design  to  be 
acceptable . 
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.  Side  constraints  -  upper  and  lower  bounds  in  a 

design  variable . 
Assuming  that  the  FORTRAN  analysis  program  has  been 
developed  and  a  particular  objective  function  has  been 
selected,  the  general  optimization  problem  can  be  stated 
as  [Ref.  20]: 

Find  the  vector  of  design  variables, X  ,  to 

Minimize     I"  (Xy 

Subject  to  the  constraints: 


(186) 


(Sjfly^O      j=1fNC0tf  (187) 

Hj(x)=0       i*i,NEQ  (188) 

VLBj^Xi^VUB;   1*1.  NOV  C189) 

where        X  =  the  vector  containing  the  set  of  independent 

design  variables. 
fCx)     -   the  objective  function  to  be  minimized. 
<£ Ax)    =  inequality  constraint  (NCOAi  is  the  number 
of  such  constraints) . 
Wj(x)  =  equality  constraint  (  NEC?  is  the  number 
of  such  constraints). 
VLB-/VU8'L=  lower  and  upper  bounds,  respectively, 
on  the  design  variables. 
If  all  inequalities  of  Eqs .  (187)  and  (189)  are  satisfied, 
the  design  is  said  to  be  feasible  if  any  constraint  is  not 
satisfied,  the  design  is  infeasible.   If  the  objective  function 
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is  a  minimum  and  the  design  is  feasible,  it  is  said  to  be 
the  optimal  design. 

In  order  to  start  the  optimization  algorithm,  the 
initial  set  of  design  variables, X  ,  must  be  specified.   It 
is  desirable,  but  not  essential,  that  the  initial  design 
variables  provide  a  feasible  solution.   The  optimization 
algorithm  will  then  proceed  in  an  iterative  fashion  using 
the  following  relationship 

4  +  1  _  4  -  <J 

where    %     =  the  iteration  number. 

<x  =  scalar  quantity  which  defines  the  move  in  the 

search  direction. 
5  =  vector  search  direction  which  will  reduce 
the  objective  function  (useable  direction) 
without  violating  constraints  (feasible 
direction) . 
To  solve  this  problem,  the  optimization  program 
COPES/CONMIN  is  used  [Ref.  18].   CONMIN  uses  the  Fletcher- 
Reeves  algorithm  for  locally  unconstrained  problems  [Ref.  20] 
and  Zoutendijk's  method  of  feasible  directions  (modified  to 
improve  efficiency  and  reliability  and  to  deal  with  designs 
which  do  not  initially  satisfy  all  the  constraints)  for 
locally  constrained  problems  [Ref.  21] . 

However,  CONMIN  does  not  handle  equality  constraints 
directly,  but  rather  by  means  of  penalty  parameters.  To 
achieve  this,  the  objective  function  is  augmented  as  follows 
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[Ref .  19]  : 

F'ffi-FlV-Kg  Hj  (190) 

and  the  equality  condition  of  Eq.  (188)  is  treated  as  an  in- 
equality constraint 

Hj(x)*0      j  =  1>'E<t> 

The  penalty  function  approach  effectively  satisfies  the 
equality  constraint  while  maintaining  the  rapid  convergence 
characteristics  of  the  CONMIN  program. 

The  numerical  optimization  problems  of  equations  (186) 
through  (190)  are  very  general,  allowing  for  any  number  of 
design  variables  and  constraints.   In  assessing  the  value  of 
optimization,  the  automated  design  provides  a  very  attractive 
approach  to  numerical  optimization;  however,  there  are  both 
advantages  and  limitations  to  these  techniques  [Ref.  20]. 

Advantages : 

.  Reduction  in  design  time. 

.  Systematic  design  procedure. 

.  Applicable  to  a  wide  variety  of  design  variables 
and  constraints. 

.  Virtually  always  yields  some  design  improvement. 

.  Not  biased  by  engineering  experience. 

.  Requires  a  minimal  amount  of  man-machine  interface. 

Limitations : 

.  Computer  times  may  increase  dramatically  as  the 
number  of  design  variables  increases.   A  practical 
limit  imposed  by  the  current  state  of  the  art  for 
most  problems  is  30  design  variables. 

121 


Optimization  techniques  have  no  stored  experience 
to  draw  upon;  the  validity  of  the  result  is 
limited  to  the  validity  of  the  analysis  program. 

The  results  of  the  optimization  are  as  correct 
as  the  analysis  program  is  theoretically  precise. 

Optimization  algorithms  used  here  cannot  deal  with 
discontinuous  functions. 

The  optimization  program  will  not  always  obtain  a 
global  design  optimum  and  may  require  restarting 
from  several  different  points  to  acquire  reason- 
able assurance  of  obtaining  the  global  optimum. 

The  analysis  program  must  be  properly  structured 
to  couple  with  the  COPES/CONMIN  optimization  code. 


C.   DESIGNATED  DESIGN  VARIABLES,  CONSTRAINTS  AND  OBJECTIVE 
FUNCTION 

To  assist  in  the  interpretation  of  the  enclosed  OTEC 
power  system  FORTRAN  analysis,  the  following  summary  identi- 
fies the  design  variables,  constraint  functions  and  objective 
function  used  in  the  analysis  and  subsequently  operated 
upon  by  the  COPES/CONMIN  optimization  code.   These  parameters 
are  all  contained  in  a  labeled  COMMON  block  in  the  computer 
code,  referred  to  here  as  "GLOBAL  COMMON."   Specific 
GLOBAL  COMMON  location  numbers  and  upper/lower  bounds  for 
operating  parameters  summarized  below  can  be  located  in 
Appendix  C. 

Design  Variables 

.  Inner  cold  pipe  diameter 

.  Inner  hot  pipe  diameter 

.  Inner  ammonia  circ  pipe  diameter 

.  Inner  ammonia  re-flux  pipe  diameter 
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Evaporator  operating  pressure 
Condenser  operating  pressure 
Outer  condenser  tube  diameter 
Outer  evaporator  tube  diameter 
Evaporator  tube  length 
Condenser  tube  length 
Condenser  tube  salt  water  velocity 
Cold  pipe  salt  water  velocity 
Evaporator  tube  salt  water  velocity 
Hot  pipe  salt  water  velocity 
Evaporator  tube  profile  pitch  ratio 
Condenser  tube  profile  pitch  ratio 

Constraint  Functions 
Operating  system  pressure  ratio 
Upper  temperature  bound  of  ammonia 
Lower  temperature  bound  of  ammonia 
Satisfactory  enthalpy  at  state  point  5 
Satisfactory  quality  at  state  point  5 
Satisfactory  condenser  tube  length 
Internal  turbine  efficiency 
Evaporator  tube  sheet  diameter 
Condenser  tube  sheet  diameter 

Objective  Function 
Cost  of  major  power  system  components 
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VIII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

1.  The  use  of  an  analysis  code  for  OTEC  power  systems 
coupled  to  COPES/CONMIN  optimization  code  provides  a  power- 
ful tool  to  design  an  optimum  power  system  for  the  desired 
net  electrical  output,  measured  against  the  objective 
function.   Such  a  design  could  permit  construction  of  higher 
capacity  systems  using  the  optimized  modules  as  substations 
of  the  total  power  plant. 

2.  The  analysis  code  coupled  to  COPES/CONMIN  provides 

an  excellent  vehicle  to  evaluate  proposed  designs  relative  to 
a  true  optimum.   Tables  1  through  4  illustrate  the  result  of 
preliminary  calculations  using  the  analysis  code  with  an 
objective  function  to  minimize  system  cost.   From  these,  the 
following  conclusions  can  be  drawn  concerning  horizontally 
oriented  aluminum  (Al-5052)  and  titanium-tubed  heat  exchanger 
power  systems : 

a.  The  cost/KW  output  is  nearly  constant  over  the 
range  of  optimum  designs  for  both  titanium  and  aluminum  tube 
heat  exchangers. 

b.  During  testing  for  feasible  plant  designs  in 
increments  of  5  MW  (net)  electrical  output,  it  was  observed 
that  a  higher  megawatt  output  plant  could  be  achieved  with 
titanium-tubed  heat  exchangers  than  for  aluminum  (Al-5052). 
For  titanium-tubed  heat  exchangers,  a  25  MW  (net)  power 
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system  is  a  feasible  design;  however,  aluminum -tubed  systems 
could  not  provide  a  feasible  design  for  the  same  output. 
Titanium  tubed  plants  failed  to  produce  a  feasible  design 
for  a  30  MW  (net)  output  power  system.  In  both  cases  of 
infeasible  design,  the  constraint  which  was  consistently 
violated  was  turbine  internal  efficiency,  set  at  90%  for 
current  state-of-the-art  design. 

c.  The  energy  conversion  and  efficiency  of  design 
of  a  turbine-generator  has  a  major  effect  on  the  overall 
system  performance  as  indicated  in  paragraph  b  above. 

d.  The  cost/KW  output  for  titanium- tubed  heat 
exchangers  is  one  third  the  cost/KW  output  for  aluminum-tubed 
heat  exchangers  using  a  30-year  life-cycle  criterion,  with  a 
10%  annual  inflation  rate  and  retubing  at  10  and  20  year 
marks  with  AL-5052  tubing. 

e.  Aluminum- tubed  heat  exchangers  have  larger  tube 
bundle  volumes,  with  volumetric  differences  between  aluminum 
and  titanium  varying  from  26.1  to  11.81  for  evaporators  and 
23.2  to  7.4%  for  condensers  over  the  range  of  net  poiver 
levels  considered.   In  both  cases  volumetric  differences 
diminish  as  the  system's  net  electrical  output  increases  to 
20  megawatts . 

f.  COPES/CONMIN  has  provided  optimum  designs  for 
each  incremental  output  power  level.   By  manipulating  the 
specified  design  variables,  subject  to  imposed  constraints, 
COPES/CONMIN  has  created  designs  whose  geometry  and  operating 
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parameters  cannot  be  scaled  on  the  basis  of  net  power  output 
(10  MW) .   Therefore,  designs  for  component  geometry  at 
increasing  power  levels  based  upon  such  simplistic  scaling 
criteria  will  not  achieve  an  optimum  design  with  respect  to 
the  cost  objective  function. 

B.   RECOMMENDATIONS 

1.  Evaluate  additional  objective  functions  including: 

a.  Minimize  heat  exchanger  volumes. 

b.  Minimize  parasitic  power  losses. 

c.  Maximize  thermodynamic  efficiency. 

d.  Maximize  net  electrical  output. 

2.  Perform  a  sensitivity  analysis  on  power  system  design 
variables  to  evaluate  their  influence  on  component  and  system 
performance.   This  allows  the  designer  to  prioritize  system 
components  which  can  provide  improvement  in  the  objective 
function  for  a  corresponding  improvement  in  component  design. 

3.  Considerable  uncertainties  are  associated  with  the 
expressions  used  to  estimate  component  performances  (two-phase 
pressure  drops,  film  coefficients,  etc.).   The  code  should  be 
tested  to  determine  the  sensitivity  of  system  design  to  these 
uncertainties . 

4.  Expand  the  code  to  include  the  use  of  enhanced  heat 
transfer  techniques  and  evaluate  the  influence  of  increased 
piping  friction  factors  on  pumping  power  requirements. 
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5.  Evaluate  proposed  OTEC  designs  using  proposed  system 
parameter  inputs,  comparing  both  the  basic  analysis  and  the 
optimization  output. 

6.  Select  other  analytical  expressions  for  heat  transfer 
coefficients  to  validate  the  performance  and  output  of  the 
existing  code. 

7.  Evaluate  the  effect  of  a  smaller  thermal  difference 
seen  by  the  power  system  and  its  influence  on  a  feasible 
design  for  a  specific  net  electrical  output. 

8.  Evaluate  the  cost  aspects  of  using  variable-pitch 
pumps  versus  fixed-blade  for  a  variable  thermal  gradient 
environment . 

9.  Evaluate  and  verify  the  influence  of  incremental 
improvements  (percent)  in  turbine  internal/adiabatic  efficiency 
with  respect  to  gross  and  net  electrical  outputs  and  compare 
with  the  results  reported  in  Ref.  16. 
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APPENDIX  A 

SAMPLE  INPUT  DATA  FOR  OTEC  ANALYSIS 

EVAPORATOR  -  HORIZONTAL 

TO  BE    O.O.  L.  000  (IN)  25.400 (MM) 

TOBE    LENGTH  tO.OOOtFT)  12.192(H) 

SW    TOBE    VEL  6.0C0(FT/S)  i.829(M/S) 

OPER    PRESSORE  1 30.000 { LBF/ 1  HZ  )  0.896(MPA) 

TOBE    MATERIAL    -    TITANIUM 

THERMAL    CGND(K)  9.500( BTU/HR.FT.F J  16. 502( W/M.C ) 

TUBE    PROFILE    -    STAGGERED    EQJI-LATERAL 

PITCH    RATIO  1.50 

ENHANCEMENT    -    PLAIN    TU3E 

CONDENSER    -    HORIZONTAL 

TUBE    0.0.  l.OOO(IN)  25.4Q01MM) 

TUBE    LENGTH  56.500(FT)  17.22KM) 

SW    TOBE    VEL  6.0Q0(FT/S)  1.829(M/S) 

OPER    PRESSORE  d9.00Q( LBF/IN2 )  0.ol4(MPA) 

TUBE    MATERIAL    -    TITANIUM 

THERMAL    CCNO(K)  9. 500 ( BTU/HR . FT.F )  16.502 (W/M.C) 

TUBE    PROFILE    -    STAGGcRED    EQOI-LATERAL 

PITCH    RATIO  1  . 50 

ENHANCEMENT    -    PLAIN    TUBE 

SALT    WATER    HOT    PIPE 

PIPE    I.D.  19.300(FT)  5.333(M) 

PIPE    LENGTH  3G0.000(FT)  ?l.440(M) 

SW    PIPE    VEL  4.500IFT/S)  i.372(M/S) 

SW    INLET    TEMP  80.000IDEG    F)  26.567(DES    C) 

SW    SALINITY  35.    0/000 

SALT    WATER    COLO    PIPE 

PIPE    I.D.  19.600(FT)  5.974(M) 

PIPE    LENGTH  3000.000{FT)  91^.400(M) 

SW    PIPE    VEL  5.500(FT/S)  1.676(M/S) 

SW    INLET    TEMP  40.000(DEG    F)  *.«*44(DEG    C) 

SW    SALINITY  35.    0/000 
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AMMONIA    CIRC    PIPE 

PIPE     I.D.  2.00CHFT)  0.610M) 

PIPE    LENGTH  L5J.00O(FT)  *5.720(M) 

AMMONIA    RE-FLUX    PIPE 

PIPE     I.D.  2. OOO (FT)  0.610 (M) 

PIPE    LENGTH  50.0C0(FT)  15.2<+0(M) 

PUMP    AND    GEN-TURB    PERFORMANCE 
EVAP    S*    PUMP 

EFFICIENCY    MECH       95.00(PCT)  4GTGR    ^8.3:(PCT) 

CONO    SW    PUMP 

EFFICIENCY    MECH       85.G01PCT)  MOTOR    98.03(?CT) 

AMMONIA    CIRC    PUMP 

EFFICIENCY    MECH       75.00(PCT)  40T0R    98.00(PCT) 

GEN-TURB    EFFICIENCIES 

GEN    MECH&ELECT  9&.60(PCT) 

TURB    MECH  99.80<PCD 

POWER    REQUIREMENTS 

NET    POWER    OUTPUT  L5.0001MW) 
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APPENDIX    B 
SAMPLE   OTEC  ANALYSIS    OPTIMIZATION   OUTPUT    DATA 

EVAPCRATCR    -    HORIZONTAL 

HT    ABSORB         2  00501 3760.01 3TU/HR )  587.60UMW) 
SW    FLOW                  334081280. 0(LBM/HR)             1 51535904. 0( KG/HR ) 

SW    TEMP    IN                 30.0G0(DEG    F)  26.567{DEG    C) 

SW    TEMP    OUT               73.724(DEG    F)  23.130(0EG    C) 

NH3    FLOW                   3788703. OCLBM/HR)  1718 519. J ( KG/HR  ) 

OPER    PRESSURE             1 30. 097( LBF /I N2 )  396.987(KPA) 

EVAP    SAT    TEMP          70.5B5{OEG    F)  21.43&OEG    C) 

OUTLET    TEMP               70.514(DEG    F)  21.397(D£G    C) 
OUTLET    QUALITY            92.00(PCT) 

NH3    PRESS    DROP               0. 162 ( LBF/ IN2 )  1.115(KPA) 
TUBE    CHARACTERISTICS 

GUTTER    DIA                     0.952UN)  24.180(MM) 

WALL    THICK                     0.025(IN)  0.639(MM) 

LENGTH                             42.182(FT)  12.357{M) 
MATERIAL    -    TITANIUM 
TUBE    PROFILE    -    STAGGERED    5QUI-LATERAL 

PITCH    RATIO  1.40 

ENHANCEMENT    -    PLAIN    TUBE 

SW    VELOCITY                  6.026(FT/S)  1.337M/S) 

T    WALL(SHELLSIDE)     7l.466(DEG    F)  21.926CDEG    C) 

FILM    TEMP                         70.990<DEG    F)  21.66KDEG    C) 

OELTA    T    3UILING            0.9521DEG    F)  9.52910EG    C) 

L.M.T.D.                                5.757(0EG    F)  3.19310EG    C) 
EVAP    EFFECTIVENESS                             0.66 T 
NR    OF    TRANSFER    UNITS                       1.093 

OVL    HT    COEF                    612 . 97 ( 8TU/HR . FT2 . F )  3430 . 58( A/ M2.C ) 

H{WATER)                     1  1I5.97(BTJ /HR.FT2.F J  5336. 71(  W/M2.C) 

H(FOULING)                3733. 54( S T J/HR. FT2 . F J  21 512 . 09(  W/M2 .C  ) 

H(METAL)                      4409.72(3TU/HR.FT2.F)  25039. 29( W/M2.C J 

H(AMMONIA)                4038.44(3TJ/HR.FT2.F )  232 1 5 . 00(  W/ AZ  .C J 
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HT    SURFACE               572405 . 00  I F T2 )  53178. 121  AZ) 

TUBE    SHEET    OIA            27.2I3(FT)  8.294(,M) 

TOT    NR    OF    TUBES  54449. 

SW    PRESS    OROP                  4.06KLBF/IN2)  2  3.  00  5  (XPA) 
MOISTURE    SEPARATOR-INSIDE    EV  AP    SHELL 

OPER    PRESSURE             1 29. 935( LBF / I N2  )  395.372<KPA) 

OUTLET    TEMP               70.333(0£G    F)  21.296(0EG    C) 

OUTLET  QUALITY    99.5G(PCT) 

NH3    PRESS    DROP          0.4 L6< LBF/ I N2 )  2.868CKPA) 

CONDENSER    -    HOKIZUNTAL 

HT    REJECT          19354723^6. 013TU/HR)  567.22WMW) 

SW    FLOW  334871552. 0<L3M/HR)  1 5 1394363. 0( KG/HR  ) 

SW    TEMP     IN                  40.000(OEG    F)  4.444(0EG    C) 

SW    TEMP    OUT               <+6.055<DEG    F)  7.308<0£G    C) 

NH3    FLCW                    3786708.  (MLBM/HR)  L  718519.  0(  KG/HR  J 

OPER    PRESSURE               o3. 1 51( LBF / I N2  )  607.78HKPA) 

CONO    SAT    TEMP          49.36KDEG    F)  9.o45(DEG    C) 

CUTLET    TEMP               49.238(DEG    F)  9.577(DEG    C) 

NH3    PRESS    DROP               0. 206( LBF/ IN2  )  l.t23(KPA) 

TUBE    CHARACTERISTICS 

OUTTER    DIA                    3.972(IN)  24.683CMM) 

WALL    THICK                    0.025UN)  0.6421MM) 

LENGTH                             57.416(FT)  L7. 500(H) 
MATERIAL    -    TITANIUM 
TUBE    PROFILE    -    STAGGERED    EQUI-LATEPAL 

PITCH    RATIO  1.40 

ENHANCEMENT    -    PL  AIM    TUBE 

SW    VELOCITY                  6.017JFT/S)  1.834M/S) 

T    WALL(SHELLSIDE)     43.33KDEG    F)  9.073OEG    C) 

FILM    TEMP                          48.785{DEG    F)  9.3251DEG    C) 

•DELTA    T    CONO                     0.907(DEG    F)  0.504(0EG    C) 

L.M.T.D.                                5.683(DEG    F)  3.157(DE3    C) 

CONO    EFFECTIVENESS  0.655 

NR    OF    TRANSFER    JNITS  1.065 
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CVL    HT    CGEF 
H{  HATER  J 

Hi  FOULING) 
H(METAL) 
H{ AMMONIA) 

HT    SURFACE 

TUBE    SHEET    DIA 


<+46.33(BTJ/HR.FT2.F  )  2537 . i 3( W/M2 .C ) 

704.35CBTU/HR.FT2.FJ  39  99.4H  W/M2.C) 

3  792.00<3TJ/HR.FT2.F  )  21 531. 73{ W/M2 .C ) 

'+393.6316TU/HR.FT2.F  )  24947, 90<  W/M2.C) 

3G53.o5(3TU/HR.FT2.F  )  17367. 63 < W/M2.C 1 

762190. 25{FT2)  70809. 69<M21 

27.194(FT)  8.2891.M) 


TOT    NR    OF    TUBES 

SW    PRESS    DROP 


52L79. 
5.636ILBF/IN2) 


3  3.36H  ,<P4) 


SALT    WATER    HOT    PIPE 

PIPE     I.D.  20.077(FT) 

PIPE    LENGTH  300.000{FT) 

SW    PIPE    VEL  4.5971FT/S) 

SW    FLOW  334081230. 0(LBM/HR) 

SW    INLET    TEMP  30.000OEG    F) 

SW    SALINITY  35.    0/000 

SW    PRESS    DROP  0.3221L3F/IN2  ) 

SALT    WATER    COLD    PIPE 

PIPE     1.0.  13.622(FT) 

PIPE    LENGTH  3000.000(FT) 

SW    PIPE    VEL  5.334(FT/S) 

SW    FLOW  334871552. 0(LbM/HR) 


SW    INLET    TEMP 
SW    SALINITY 
SW    PRESS    OROP 
AMMONIA    CIRC    PIPE 
PIPE    I.D. 
PIPE    LENGTH 
NH3    FLOW 


40.00CMDEG    F) 

35.    0/000 
0.508(LBF/IN2 ) 

2.00KFT) 
150.0G0(FT) 
3733708. 0(LBM/HR) 


NH3    PRESS    DROP  15. 0 33( LSF / IN2 ) 

AMMONIA    RE-FLUX    CIRC    PIPE 

PIPE     I.D.  2.000(FT) 

PIPE    LENGTH  50.0C0(FT) 

NH3    FLGw  1136612. OtLBM/HR) 

NH3    PRESS    DROP  9. 8 74( LBF/ IN2  ) 


6.120CMJ 
9L440(MJ 
J..40U  M/S) 
151535904.  Q<  KG/HR ) 
26.6670EG    C) 

2.217(KP4) 

5.676(M) 
914.400C Ml 
L626CM/S) 

151394363.  0(  KG/HR  ) 
4.*44(DEG    C) 

3.50KKPA) 

Q.610(M) 
45.72  J(  M) 
1713519.  0(  KG/HR  ) 
L03.6501KPA) 

0.610( M) 
15.240(  A) 
515555. 3(KG/HR) 
68.079< KPA) 
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PUMP  AND  GEN- TURB  PERFORMANCE 

EVAP  SW  PUMP 

HEAD  PRESS        9.8  98(FT) 
CAPACITY       653259. 6<GAL/MIN) 
EFFICIENCY  MECH   85.00(PCT) 

COND  Sw  PUMP 

HEAD  PRESS       2Q.753(FT) 
CAPACITY       652119. 2(GAL/MIN) 
EFFICIENCY  MECH   35.00(PCT) 

AMMONIA  CIRC  PUMP 

HEAD  PRESS      2I2.327(FT) 
CAPACITY        121  )1.  f|  GAL /.MINI 
EFFICIENCY  MECH   75.00<PCT) 

AMMONIA  RE-FLUX  PUMP 

HEAO  PRESS       33. 016(FT) 
CAPACITY         3732.4(GAL/MIN) 
EFFICIENCY  MECH   75.00{PCT) 

GEN- TURB  EFFICIENCIES 

GEN  MECHdELECT     96.60IPCT) 
TURB  MECH  99.30(PCT) 

TURB  INTERNAL      39.33(PCT) 

TURB  OUTLET  QUALITY  96.77(PCT) 


3  .  J  i  7  (  M  ) 
2472  567.  0(  LI  T/MIN) 
MOTOR    98. OCX  PCT) 

9.22  7(M) 
2468271.  J(  LI T/MIN) 
MOTOR    98.03( PCT) 

64.717(  M) 
45805. 0( LI T/MIN) 
MOTOR    99.00( PCT) 

1 1 .  5  3  7  (  M ) 
14127.  II  LIT/MIN) 
MOTOR    98.  3  0 (PCT) 


POWER  REQUIREMENTS 

TUR8-GEN  GrfOSS  276b3 .3 13 ( HP ) 

EFFICIENCY  LOSSES 
EVAP  SW  PUMP     19o4.851(HP) 
COND  SW  PUMP     4129.313(HP) 
NH3  CIRC  PUMP     541.714IHP) 
NH3  RE-FLUX  PUMP   29.097(HP) 


20.6  33 (MW) 
0  .  5  5  9  (  MW  ) 
1.496(MW) 
3.143(MW) 
0.412 (MW) 
0.022 (MW) 


NET    POWER    OUTPUT 


15.000 (MW) 


PERCENT    PARASITIC    POWER 
THEMODYNAMIC    CYCLE    EFFICIENCY 


2 4. 59 (PCT) 
2.65( PCT) 
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COST    QF    COMPONENTS 
EVAPORATOR 
CONDENSER 
GEN-TUKBINE 
GENERATOR 
EVAP    SW    PUMP 
COND    SW    PUMP 
NH3    CIRC    PJMP 
NH3    RE-FLUX    PUMP 


8223672.  00(  DOLLARS  ) 
8667154. 00( DOLLARS) 
L578776-00CD0LLARSJ 

937205. 06(D0LLARS) 
653332. 94(D0LLARS) 
652298.06( DOLLARS) 
1 36824. 94( DOLLARS) 
64443. 34( DOLLARS) 


CPTIMUM    COST 
COST    PER    NET    Kw    OUTPUT 


20349216. 00( DOLLARS) 
L339.951 DOLLARS) 
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APPENDIX    C 
SAMPLE    COPES    OPTIMIZATION  AND   SENSITIVITY   ANALYSIS    DATA 

$BLQCK    A    (TITLE    CARD) 

CCEAN    THERMAL    ENERGY    CONVERSION     (OTEC)     POWER    SYSTEM 

SBLQCK    3     (PROGRAM    CONTROL    PARAMETERS) 

2,16,16 

2         16         16 
S6L0CK  C  (INTEGER  OPT  CONTROL  PARAMETERS) 
5,0,0,5 

5  0  0  5 

$BLGCK  0  (FLOATING  PT  OPT  PROG  PARAMETERS) 
0-0 

0.0 
0.0 

0.0 
$6L0CK    E     (TOT    NR    DESIGN    VAR, DESIGN    03J'  IDENT    AND    SIGN) 
16,27,-1.0 

16  27  -1.0 

JBLOCK    F     (DESIGN    VARIABLE    BOUNDS, INIT    VALJES    u    SCALE    FACTOR) 
1.0,1.0  +  20 


1.0 

1.0+20 

1.0,1.0+20 

1.0 

1.0+20 

1.0,1.0+20 

1.0 

1.0+20 

1.0,  1.0  +  20 

1.0 

1.0+20- 

85.0,1*3.0 

35.0 

148.0 

85.0,148.0 

35.0 

148.0 

0,5,2.5 

0.5 

2.5 

0.5,2.5 

0.5 

2.5 

10.0,1.0+20 

10.0 

1. G+20 

10.0,1.0+20 

10.0 

1.0+20 

2. Ot  10.0 

2.0 

10.0 

2.0,  10.0 

2.0 

10.0 

2.0,  10.0 

2.0 

10.0 

2.0,10.0 

2.0 

10.0 

1.4,3.0 

1.4 

3.0 

1.4, 3.0 

1.4 

3.0 

16L0CK    G    (DES1 

;SN    VARIABLE 

IDEN 

1, 1,  1.0 

1 

1 

1.0 

2,1,  1.0 

2 

2 

1.0 

3,3,1.0 

3 

3 

1.0 

4,4,  1.0 

4 

4 

1.0 

5,5,  1.0 

5 

5 

1.0 
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6,6, L.O 
7,7,  1.0 
8,3,  1.0 
5,9,  1.0 


5 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


1.0 
1.0 
L.O 
L.O 
1.0 
L.O 
1.0 
1.0 
L.O 
L.O 
1.0 


10,10,1.0 

10 
11,11,1.0 

11 

12,12,1.0 

12 
13,  13,1.0 

13 
14, 14,1.0 

14 
15,15,1.0 

15 
16,16,1.0 

16 
SBLOCK    H    INR    Of    CONSTRAINED    PARAMETERS) 
5 

5 
SBLQCK     I     (CONSTRAINT    I  DENT    AND    BOUNDS) 
17 

17 
l.,0.0, 3. ,0.  0 

1.  0.0  3.  0.0 

18,20 

13  20 

0.1,0.0,1.0+20,0.0 

0.1  0.0  i.0+20  0.0 

21,23 

21  23 

0.,0.0,  1.0+20,0.0 

0.  0.0  1.0+20  0.0 


24 


24 


30. ,0.0,90. , 0.0 

30.  0.0  90.  0.0 

25,26 

25  26 

10. ,0.0* 1.0+20,0.0 

10.  0.0  1.0+2  0  0.0 

S8L0CK    P    (SENSITIVITY    OBJECTIVES) 
27,0 

27  0 

1,2, 3, 4, 5, 6, 7, 3, 9, 10, 11,  12, 13, 14, 15, 16, 17, 13, 19, 20, 


12  3  4 

9  10  11  12 

17  18  19  20 

25  26  27 

21,22,23,24,25,26,27 

6  7  3 

14  15  16 

22  23  24 

$BL0CK    Q    (SENSITIVITY    VARIABLE    BOUiNOS) 
1,6 

1  6 
18., 10. ,15., 13., 20. ,25. 

13.  10.  15.  18. 

2,5 

2  6 
18. ,10. ,15. ,18.  ,20.  ,25. 

13.  10.  15.  13. 

3,6 


5 
13 
21 


20. 


ZO. 


25. 


25. 


158 


3  o 

2.  • 0«5t 1 . » 2.  f  3.  , 4. 

2.        0.5         L. 
4,6 

4-        ^   6 

2.        0.5         1. 
5,6 

5  6 

130.  ,12 3., 12 9., 130. ,131. ,132. 

130.       123.        129. 
6,6 

6  6 
89.,  67.  ,33.  ,t)9.  ,  90.  ,91  . 

39.  37.        83. 
7,6 

7  6 
1.  ,  .5, . 75,1. , 1.25,  1.5 

1.         .5        .75 
8,6 

8  6 
1. ,.5, .75,1. ,1.25,1.5 

1.         .5        .75 
9,6 

9  6 
40.,  30. ,3  5. ,40., 45., 5  0. 

40.  30.        35. 
10,6 

10  6 
55. ,45. ,50. i55.f60.f65. 

55.        45.        50. 
11,6 

11  6 
7.,5.»6.,7.,8.,9. 

7.         5.         6. 
12, o 

12  6 
4.,2.,3.f4.f5.f6. 

-+.         2.         3. 
13,6 

13  6 
7«f5.f6. »7.»8.»9. 

7.         5.         6. 
14,6 

14  6 
4. ,2. ,3. ,4., 5. ,6. 

4.         2.         3. 
15,6 

15  6 
l.o, 1.4, 1.5, 1.6, 1.7, 1.8 

1.6        1.4        1.5 
16,6 

16  6 
1.6, 1.4,1.5,  1.6,  1.7,1. 3 

1.6        1.4        1.5 
END 


2. 


2. 


13  0. 


89. 


1. 


i. 


40. 


5d. 


7. 


7. 


1.6 


1.6 


3. 


3. 


131. 


90. 


1.25 


1.25 


45. 


60. 


3. 


5. 


3. 


3  . 


1.7 


1.7 


4. 


4. 


132. 


91. 


1.5 


1.5 


50. 


65. 


9. 


6. 


9. 


6. 


1.8 


1.3 
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